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VIRTUAL PROTOTYPING AND TESTING FOR 
MEDICAL DEVICE DEVELOPMENT 

BACKGROUND 

5 

1 . Field of the Invention 

This invention relates to systems and methods of developing better-designed medical 
devices, specifically, intracorporeal medical devices and particularly cardiovascular stents and 
endovascular grafts. 

10 2. Background and Description of Related Art 

Atherosclerotic vascular disease is a significant health problem facing the world 
population today. Atherosclerosis results in two primary types of lesions— occlusive and 
aneurysmal, with the aorta being the primary site of aneurysmal disease. Occlusive disease is a 
process in which a vessel lumen becomes narrowed and the blood flow restricted. Occlusive 

1 5 disease is typically associated with plaque buildup on the vessel wall or a biological response to 
vessel injury. One approach to treatment of occlusive disease involves placing a stent inside the 
vessel to act as a structural scaffold and hold open the vessel, and also possibly to provide local 
drug delivery or local radiation treatment. Aneurysmal disease is a process in which a vessel 
dilates under the influence of hemodynamic pressure, and may ultimately lead to rupture of the 

20 vessel and severe internal bleeding. One approach to treatment of aneurysmal disease involves 
placing a TPEG (transluminally placed endovascular graft, or "stent graft") across the aneurysm, 
excluding the aneurysm from hemodynamic pressure and thereby reducing or eliminating the risk 
of rupture. Examples of such grafts can be found in co-pending U.S. Patent Application Serial 
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No. 09/133,978, filed August 14, 1998 by Chobotov, which is hereby incorporated by reference 
herein in its entirety. 

A TPEG is an endovascular prosthetic device that lines the interior of an artery to provide 
flow path integrity and structural support to the damaged or diseased blood vessel. TPEGs are 
5 sometimes called "stent grafts" because they were originally created using combinations of stents 
and synthetic vascular graft segments. TPEGs are delivered to a blood vessel location in a 
compressed state, through an incision, and are then deployed at the location of concern. 

The current development process of TPEGs and medical devices generally, usually 
involves the reiterative and sequential steps of designing, fabricating the prototype, and testing 

1 0 the prototype until the required performance specifications are met. Fabrication of the prototype 
entails the building of the actual medical device, e.g., a TPEG. Testing can involve animal 
testings, human clinical trials, stress, strain, and deformation testing, and the like. Stents, TPEGs 
and other medical devices have suffered from long development times and from design 
deficiencies discovered late in the development and testing process. Thus, the development of 

1 5 improved medical devices could be significantly accelerated if design deficiencies could be 
identified earlier, before committing to lengthy laboratory testing, animal studies, and human 
clinical trials. A system that enables early evaluation of many aspects of device performance in 
vivo, and is applicable to development of stents for occlusive disease, TPEGs for aneurysmal 
disease, and other medical devices is highly desirable. 

20 In designing a TPEG, several factors must be taken into account, such as the structural 

integrity of the TPEG, the prevention of perigraft leaks, the need for a more easily-controlled 
TPEG deployment to allow a more precise positioning of the TPEG, the kink resistance of the 
TPEG, the morphology of the arterial walls, the relatively large size and lack of TPEG flexibility 
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in the undeployed configuration (which can create difficulties in passing the TPEG from its 
insertion site to its deployment site), and the like. In vivo boundary conditions and forces, 
particularly dynamic or static cyclic in vivo forces, and the material properties of a TPEG are 
also important factors. Taking these factors into consideration during virtual testing and 
5 development of a medical device generates a more accurate assessment of the maximum stresses, 
strains, and deformations, over time that may potentially be handled by a medical device such as 



In designing a stent, several factors must be considered including radial force, crush 
resistance, flexibility (in both the compressed and the deployed configurations), fatigue life, and 

10 tissue intrusion through open stent cells. A system that allows rapid evaluation of these and 
other characteristics of a stent design before hardware prototypes are constructed, thereby 
reducing the cost and time required for development and also expanding the designer's capability 
to explore more exotic designs and possibly discover new and more advantageous stent designs 
within a given budget and timeframe is highly desirable. 

1 5 Thus, systems and methods which allow accurate virtual testing of a medical device 

design with respect to one or more of the above noted factors, in addition to other factors not 
specifically enumerated, without the need for an actual prototype of the design, are needed. Such 
systems and methods can reduce the cost of medical device development and increase the safety 
and efficacy of the designs. 



The invention provides a system and method for developing better-designed medical 
devices and particularly cardiovascular stents and endovascular grafts. The system comprises a 
Geometry Generator, a Mesh Generator, a Stress/Strain/Deformation Analyzer, and, optionally, a 



a TPEG. 



20 
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Visualization tool. The invention may obtain anatomic data from 3D volumetric data. In other 
embodiments, the invention utilizes an idealized anatomical feature, an in vitro model, or no 
anatomical feature at all. 

In one embodiment, the Geometry Generator receives three-dimensional volumetric data 
5 of an anatomical feature and accordingly extracts the surface points of such data, which in turn is 
received by the Mesh Generator. In another embodiment, the Geometry Generator based on 
algorithms available in such Geometry Generator software generates an output that is directly 
received by the Mesh Generator. Using the output generated by the Geometry generator and the 
geometric model of a candidate medical device, the Mesh Generator generates a mesh or a finite 

10 element model incorporating either the anatomical feature or in vitro model and candidate 
medical device. In an embodiment where no anatomical feature is used, a mesh only 
incorporating the candidate medical device is generated. The Stress/Strain/ Deformation 
Analyzer then receives the mesh and the material models, the loads and/or displacements placed 
on the anatomical feature or in vitro model, if applicable, and the candidate medical device. 

1 5 Using stress and strain deformation analysis, particularly non-linear analysis, the 

Stress/Strain/Deformation Analyzer simulates and analyzes the potential in vivo stresses, strains, 
and deformations or motions of the candidate medical device. Such strains, stresses, and 
deformations may optionally be displayed using a Visualization tool. 

Various embodiments of the invention can be used to provide a variety of useful 

20 functions and capabilities to those who design, manufacture and use medical devices. 

Specifically, embodiments of the invention may be used to model anatomical features or 
anatomical environments dynamically. As a result, a computer generated model of a medical 
device, or the like, may be virtually placed or deployed within the anatomical model to measure 
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the response of the device to the environment. The dynamics of the computer generated model 
of the anatomical features or environment can be accelerated dramatically such that large 
numbers of normal biological cycle, such as a heartbeat, can be imposed upon the computer 
generated medical device model in a relatively short period of time. 
5 This gives medical device designers the ability to virtually test a proposed design in a 

short period of time relative to the time it would take for a similar number of dynamic biological 
cycles in vivo. Thus, the iterative process of device design and testing of designs is accelerated 
and improvements in medical device technology can be achieved at a quicker rate. Further, 
embodiments of the invention can be used to vary and test material properties of medical device 
1 0 components over a broad range in a short period of time using the non-linear modeling 

cababilities of the embodiments. This capability can be used to select materials having optimal 
properties for producing the safest and most efficacious designs within a given set of design 
parameters. 

Another benefit of embodiments of the invention is directed to varying material and 
1 5 configuration properties of models of anatomical features such that a simulation of testing of a 
given device could be performed in a large number of patients, as might be carried out in a large 
scale clinical trial. If the statistical variation of tissue parameters of a given anatomical feature is 
known for a given patient population, a medical device model could be tested in anatomical 
models which vary over such a given range. In this way, a large scale clinical trial could be 
20 modeled with embodiments of the invention, at least as to certain performance parameters, 

without the need for large numbers of actual patients being subjected to clinical testing. The data 
generated from such a clinical trial modeling exercise could be used to produce or refine the 
design of a medical device such that it performs optimally over a broad range of anatomical 
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environments. The design could be refined using such data to improve robustness and 
adaptability of the medical device design. 

Also, it is possible to use embodiments of the invention to identify failure modes of given 
medical device designs when such designs are subjected to dynamic mechanical and chemical 
5 forces. By identifying the cause of failure in a design, the "weak link" in the design can be 
pinpointed and necessary corrections to materials or configuration made in order to obviate the 
problem. It is also possible to test theories of failure experienced during in vivo clinical testing 
using embodiments of the invention. In other words, if an in vivo clinical failure of a medical 
device should occure, there may be one or more theories postulated as to the cause of the failure, 

10 particularly in a situation where multiple components of a device have failed and it is not clear 
from the clinical data which failure ocurred first, or if an initial failure of one component of the 
device precipitated subsequent failure of other components of the device. The dynamic modeling 
capabilities of embodiments of the invention can allow rapid testing of multiple theories as to the 
timing and causation of complex failure modes and quickly determine which of the postulated 

1 5 theories is correct. 

In addition, the dynamic, non-linear analysis modeling capabilities of embodiments of the 
invention allow a physician, who is responsible for use or implementation of a medical device, to 
more accurately choose a proper size or type of medical device based on a specific patient's 
anatomy. Such is the case when a specific patient's anatomy or anatomical feature is 

20 substantially duplicated by a computer model of an embodiment of the invention generated from 
3-D volumetric image data, or the like. A large number of sizes or types of virtual medical 
devices can then be placed and tested within the patient's specific anatomical feature to 
determine optimum safety and efficacy of the design choice. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 illustrates a block diagram representation of a virtual prototyping system having 
features of the present invention. 

Figure 2 illustrates a block diagram showing data received by an embodiment of a 
5 Geometry Generator and a Mesh Generator in accordance with the present invention. 

Figure 3 illustrates a block diagram representation of another embodiment of a system of 
the present invention. Figure 4 illustrates a block diagram showing data received by a 
Stress/Strain/ Deformation Analyzer. 

Figures 5 A-5M contain an exemplary text of a command file that is read by a Mesh 
10 Generator, such as TRUEGRID, to conduct a component-level analysis of a stent, without the 
option for simulating deployment into CT-based anatomy. 

Figures 6A-6F contain an exemplary text of a command file read by TRUEGRID for a 
simulated TPEG graft deployment in a proximal aortic neck to generate a mesh incorporating 
both an anatomical feature and medical device and to output files that are read by a 
1 5 Stress/Strain/Deformation Analyzer. 

Figures 7A-7C contain an exemplary include file used by the command file listed in 
Figures 6A-6F. 

Figures 8A-8L contain another exemplary command file read by TRUEGRID used in the 
virtual prototyping system of the present invention for simulating stent deployment into an 
20 anatomy from CT data, as opposed to a stent graft deployment. 

Figures 9A and 9B illustrate a process to develop better-designed medical devices, 
particularly TPEGs, in accordance with an embodiment of the present invention using 3D 
volumetric data. 
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Figure 10 illustrates a process to develop better-designed medical devices using in vitro 
anatomical features. 

Figure 1 1 illustrates the use of an embodiment of the present invention as a physician 
preprocedure planning tool. 
5 Figure 12 contains a representation of one simulation display of a cutaway lateral view of 

a vascular stent in the infrarenal aorta just proximal to an abdominal aneurysm. 

Figure 13 is a block diagram representation of one of the computers illustrated in Figure 

1. 

DETAILED DESCRIPTION 

10 The following detailed description illustrates an embodiment of the invention by way of 

example, not by way of limitation of the principles of the invention. Various embodiments of the 
invention will be described by way of illustration with reference to various software tools, but it 
should be understood that other software tools that have a comparable capabilities of the 
mentioned tools may be used and other medical device aside from TPEGs may also be developed 

1 5 using this invention. In addition, although the invention is discussed in the context of prosthesis 
and specifically endovascular grafts, this is in no way meant to limit the scope of the invention. 

Systems and methods of embodiments of the invention are suitable for the development 
and testing of medical devices including those for therapeutic, diagnostic, monitoring and the 
like purposes. In general, any device that interacts inside a patient's body may be better 

20 developed and tested with the systems and methods of embodiments of the present invention. 

Embodiments of the present invention are also well suited for development and testing of 
intracorporeal devices or prosthesis that generally have an acute interaction with anatomical 
features of a patent, A list of such devices, which is in no way exhaustive, could include 
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endovascular grafts, stents, pacemakers, artificial joints, artificial tendons, heart valves, artificial 
limbs, orthopedic hardware, surgical equipment such as sutures, staples, etc., and the like. 

Embodiments of the present inventions are particulary well suited for the development 
and testing of devices for use in the vascular system or other bodily systems that have stresses, 
5 strains, and deformations which are dynamic, or quasi-static, and cyclic in nature, e.g., the 
rhythmic pulsing of the arterial system resulting from variations in blood pressure from the 
patent's beating heart and the resulting cyclic dynamic or quasi-static stresses, strains, and 
deformations these variations impart on the patient's arteries and medical devices disposed 
therein or thereon. 

10 Embodiments of the present invention are also suitable for development and testing of 

interventional medical devices, which have only transient or temporary contact with the 
anatomical features of a patient. Illustrative examples of such devices can include catheters, 
balloons, atherectomy devices, guidewires, and the like. 

Figure 1 is a block diagram showing one embodiment of a virtual prototyping system 105 

1 5 for analyzing the use of a medical device constructed in accordance with an embodiment of the 
present invention. Figure 1 shows that a Geometry Generator 120 receives CT scan or MRI Data 
1 10 as input. The Geometry Generator 120 then processes the CT scan or MRI data and outputs 
data, which are then received by the Mesh Generator 130 as input. The Mesh Generator, in 
addition to receiving the output of the Geometry Generator 120, also receives a Medical Device 

20 Model data 140 as input. The Medical Device Model 140 contains the geometry (geometric 
shape or geometric model) of the candidate medical device. Such model may be the complete 
candidate, a portion, or an element of the candidate medical device. Similarly, a portion or an 
element of the anatomical features, not the entire anatomy scanned, may be received by the Mesh 
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Generator 130. The Medical Device Model may be created by a computer-aided-design (CAD) 
software application and stored as a CAD data file. Examples of suitable CAD software 
packages include I-DEAS (available from SDRC, Inc. of Milford, Ohio) and CATIA (available 
from International Business Machines Corporation), however, any other suitable application 
5 could be used. The Medical Device Model could also, for example, be created through contact or 
non-contact three dimensional measurement/imaging of a physical device or model. In another 
embodiment, the medical device model 140 is created within the Mesh Generator 130 module 
itself. 

In addition, although the embodiment of Figure 1 contemplates the use of CT or MRI 
10 volumetric data 1 10 as input, volumetric input could also be generated from any other suitable 
source, including other imaging system sources such as ultrasound imaging systems, beta scan 
imaging, radionuclide scanning, thermography and the like. Anatomical volumetric input data 
could also be artificially fabricated from idealized versions of anatomical features, which may be 
initially obtained from CT-data and modified, or be created manually by modeling such idealized 
1 5 version. These could be created to test medical devices within anatomical features having 

specified characteristics. For example, it may be desirable to test a medical device in an aorta 
having two distended sections caused by aortic aneurysms, which are separated by a non- 
distended portion of the aorta. Input data representing such an anatomical feature could be 
generated by manually entering data known to wholly represent such an anatomical feature. 
20 Alternatively, input data representing such an anatomical feature could be constructed by 

manually entering data corresponding to portions of CT, MRI or other imaging created data of 
actual patient aortas. 
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The output of the Mesh Generator 130 is then received by the Stress/Strain/Deformation 
Analyzer 160. The Stress/Strain/Deformation Analyzer 160 also receives Materials Model data 
170 and Load data 150 as input, which may also be outputs of the Mesh Generator 130. The 
output of the Stress/Strain/Deformation Analyzer 160 comprises the medical device performance 
5 data for evaluation, which may then be received by the Visualization tool 180 as input. The 

Visualization tool 180 in turn displays, through animation or visual representations, the predicted 
stresses, strains, and deformations on the candidate prosthesis 'Virtually in vivo." 

In an embodiment of the invention, the Geometry Generator 120 is a custom-developed 
software tool or the MIMICS software from Materialise NV (with offices in Ann Arbor, 
10 Michigan, USA); the Mesh Generator 130 is TRUEGRID® of XYZ Scientific Applications, Inc. 
(Livermore, California, USA); the Stress/Strain/ Deformation Analyzer 160 is a modified version 
of NIKE3D or DYNA3D available from Lawrence Livermore National Laboratory (LLNL); and 
the Visualization tool 180 is the GRIZ visualization software, also developed by LLNL. 

The unique combination of tools, data, and processing techniques as described herein in 
1 5 conjunction with the preferred embodiment provides a more accurate in vitro representation of 
anticipated in vivo forces exerted on medical devices and thereby reduces cost and time in the 
fabrication and testing of prototypes. 

The various systems or components 120, 130, 160, 180, inputs (e.g., via files), and 
outputs (e.g., via files) of the present invention may be contained in one or in a plurality of 
20 computers. Thus, the Geometry Generator may be contained in one computer, while the 

Stress/Strain/Deformation Analyzer and the Visualization tool are run and contained in a separate 
computer. Furthermore, the inputs need not directly be received by the receiving system, e.g., 
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through a network transmission. The outputs for example, of the Geometry Generator may be 
stored in a floppy disk and read by a Mesh Generator via that floppy disk. 

Figure 2 shows the data flow for an embodiment of a Geometry Generator 120 of Figure 
1 in detail. The Geometry Generator 120 receives as input the CT scan, MRI data, or other three- 
5 dimensional (3D) volumetric data 110. It is preferred that data from CT scans or MRIs be used 
in this invention because they provide a 3D volumetric representation of patient anatomy and 
blood vessel morphology, including complex atherosclerotic plaque distribution within the flow 
lumen. This type of data thus provides an accurate representation, for example, of the 
environment on which a medical device, for example, a TPEG will be placed. The CT and MRI 

10 equipment that is used to capture such 3D volumetric data are those that are readily available. 

Certain researchers and scientists in the biological sciences have at their disposal a wealth 
of voxel data. A voxel is the unit of CT or MRI reconstructions, represented as a pixel in the 
display of the CT scan or MRI. Well-established methods to extract triangular surface 
representations (hereinafter referred to "surface points") from these voxel data using criteria such 

15 as variation in density are available. An embodiment of the Geometry Generator 120 first 
extracts the surface points, at step 220, from the CT scan or the MRI image data (e.g., 
segmentation, contour based, or 3D approach). A CAD software is then used to generate the 
Geometric Model 230 of the anatomy scanned using the extracted surface points. The extraction 
of surface points can be implemented by writing a software program that implements the 

20 techniques stated above or by available software programs. An example of a software program 
that generates surface points based on CT scan or MRI data is PREVIEW from Medical Media 
Systems. 
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The output of the Geometry Generator 120 is in the form of an Anatomy Model 240, 
which contains the geometric model of the anatomy scanned. The Anatomy Model 240 and the 
Medical Device Model 140 (containing the geometric model of the candidate medical device) are 
then received by the Mesh Generator 130 as input (usually as CAD files). The anatomy model 
5 may be a portion or an element of the anatomy scanned. Similarly, the medical device model 
may be a portion or of an element of the candidate medical device. This is useful for analyzing 
the interaction between a portion of a candidate device, such as a proximal stent in a TPEG, and 
a certain anatomical feature, such as tissue. The Mesh Generator 130 then generates a finite 
element model incorporating both the anatomy model, whether idealized or actual, and the 

1 0 medical device model as represented by box 250. 

In one embodiment, the geometric models of the anatomy and the medical device are 
created using CAD software. Generally, the geometric models are stored in the Initial Graphics 
Exchange Specification (IGES) format that is an industry-standard graphic file format for CAD 
systems. Because of its wide-use, many FEA software tools read and utilize the IGES format. In 

1 5 another embodiment, the geometric models are created directly in the Mesh Generator. 

The Mesh Generator 130 in accordance with an embodiment of the invention is 
TRUEGRID®. TRUEGRID is a 3-D finite modeling and analysis tool that generates meshes or 
finite element models. It is a software that tessellates a geometric model into hexahedron brick 
elements and quadrilateral shell elements, creating a mesh or a grid. A FEA mesh generating 

20 tool, such as TRUEGRID, uses the anatomy model 240 and medical device model 140 created by 
a CAD software to generate a mesh. In another embodiment of a Geometry Generator 120 (not 
shown in the figures), the Geometry Generator is a software tool that interfaces between scanner 
data, such as CT, MRI, and technical scanner data, and Rapid Prototyping, CAD, or Finite 
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Element analysis data. Such software tools typically generate surface points from such scanner 
data, which are then converted into STL (stereolithography), slice files, and/or IGES files, which 
may then be read by the Mesh Generator 130 as input. An example of such a Geometry 
Generator 120 is the "Materialise Interactive Medical Image Control System" (MIMICS) 
5 available from Materialise, referred to above. The output of the MIMICS program, for example, 
may be directly read and processed by the Mesh Generator 1 30. Thus, steps 220 and 230, 
illustrated in Figure 2, are not necessarily implemented by this alternative embodiment of the 
Geometry Generator 120. 

Figure 3 is a block diagram showing another embodiment of a virtual prototyping system 

10 105. Figure 3 is similar to Figure 1, except that the anatomical feature is not obtained from a 3D 
volumetric data, such as a CT scan. Rather, an in vitro model of the anatomical feature is 
presented for analysis. For example, instead of a CT-scan artery, the system analyzes the 
stresses, strains, and deformations of a medical device deployed in a latex tube, which represents 
the artery or the in vitro model. Such in vitro model may be a CAD file that is read by the Mesh 

1 5 Generator 1 30 or in another embodiment created within the Mesh Generator itself. Alternatively, 
an idealized anatomical feature may be created through this embodiment. In another 
embodiment of the invention, not shown in the figure, the system may do a component or 
element analysis of a proposed medical device, without the incorporation of either an anatomical 
feature or in vitro model. 

20 Figure 4 is a block diagram showing in detail the data flow of the Stress/Strain/ 

Deformation Analyzer 160, which preferably is a non-linear finite element modeling software 
application such as DYNA3D or NIKE3D. The Stress/Strain/Deformation Analyzer receives a 
mesh incorporating both the medical device and the anatomy scanned (idealized or actual), a 
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mesh incorporating both the medical device and in vitro model, or a mesh incorporating just the 
medical device model 250. A portion of the medical device, in vitro model, or the anatomy 
scanned may be used. The Stress/Strain/Deformation Analyzer 160 also receives the Materials 
Model 170, and the Load 150 on the applicable structures (e.g., TPEG and artery or just on the 
5 medical device) to generate an output used by the Visualization tool 180. In the preferred 
embodiment, the Materials Model 170 and the Load 150 are read by TRUEGRED through a 
command file (further discussed below). Thus, the outputs of TRUEGRID (the Mesh Generator) 
do not only include the finite element model 250 of the mesh incorporating both medical device 
and anatomy scanned, mesh incorporating both medical device and in vitro model, or a mesh 

10 containing only the medical device, but the materials model 170 parameters as well as load 150 
information. This reduces the number of code changes, if necessary, within DYNA3D or 
NIKE3D, or the manual entry of input values to be read by DYNA3D or NIKE3D. 

DYNA3D is a general-purpose, explicit, three dimensional, finite element program for 
analyzing and simulating the large deformation dynamic response of inelastic solids and 

1 5 structures. DYNA3D and NIKE3D implement a number of material models, for example, 

including elastic, orthotropic elastic, and kinematics/isotropic plasticity. NIKE3D is a general- 
purpose nonlinear implicit, three-dimensional, finite element program for analyzing and 
simulating the finite strain and static and dynamic response of inelastic solids, shells, and beams. 
FEA Stress/Strain/Deformation Analyzers, such as DYNA3D and NIKE3D, are capable 

20 of analyzing and simulating sliding interfaces, body force loads due to base acceleration, body 
force loads due to spinning (geometry-dependent), concentrated nodal loads, pressure boundary 
conditions (geometry-dependent), and displacement boundary conditions. 
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The Materials Model 170 is the numerical representation of the material characteristics of 
the medical device, the anatomy, and/or the in vitro model being analyzed- Loads include 
pressures, displacement, forces, and deformations. Using the mesh 250, the Materials Model 
170, and the Load 150, the Stress/Strain/Deformation Analyzer 160 then analyzes and simulates 
5 the non-linear stress, strain, and deformation over time such as on a medical device (e.g., a TPEG 
and the arterial wall). The Stress/Strain/Deformation Analyzer in accordance with an 
embodiment of the present invention utilizes non-linear analysis (e.g., using non-linear formulas) 
or linear analysis to simulate and to analyze the non-linear static or dynamic behavior in the 
structure. 

10 In Figure 4, the Materials Model 170 is directly received by the Stress/Strain/ 

Deformation Analyzer 160. Another way to have the materials model be received by the 
Stress/Strain/Deformation Analyzer 160 is by modifying the source code of DYNA3D and 
NIKE3D, e.g., by hard-coding the materials model into the source code itself. Similarly, if the 
source code of the geometry generator, Mesh Generator, Stress/Strain/Deformation Analyzer, 

1 5 and/or Visualization tool are available, inputs as shown may be incorporated, for example, by 
actually hard-coding the input parameters into the source code or by changing certain equations 
in the code itself. 

Once the Stress/Strain/Deformation Analyzer 160 has analyzed the stresses, strains, and 
deformations on the medical device, the Visualization module 180 (in Figure 1) can then receive 
20 the output of the Stress/Strain/Deformation Analyzer to visually display the resulting stresses, 
strains, and deformations 190. 

Generally, the numerical output of the Stress/Strain/Deformation Analyzer 160 may also 
be analyzed to determine the stresses, strains, deformations on the medical device without using 
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the Visualization tool 180. Using the Visualization tool 180, however, facilitates such 
determination because the stresses, strains, and deformations are shown via a graphical and 
visual display. A virtual prototyping or simulation of a medical device design, rather than plain 
numerical output data, is thus produced. 
5 In an embodiment, the Visualization tool 180 is provided by the above-referenced GRIZ 

software application. GRIZ is an interactive software for visualizing FEA results on three- 
dimensional unstructured grids, and calculates and displays derived variables from FEA software 
tools such as DYNA3D, NIKE3D, and TOPAZ3D (also developed by LLNL). GRIZ provides 
display control of the mesh materials on an individual basis, allowing the user to concentrate on 

10 the analysis and visually focus on important subsets of the mesh, and provides the ability to 
animate the representation over time. 

GRIZ uses the Silicon Graphic Inc. (SGI) Graphics Library (GL) or Open GL for 
rendering and the "Motif widget" toolkit for its user interface. In order to compile and run GRIZ, 
both of these libraries are required. GRIZ can be used on SGI workstations as well as on SUN 

1 5 and other workstations using commercial GL emulation software. 

Considering the visual result on the screen display 190, a user may then compare the 
candidate medical device as designed against selected performance requirements. If the selected 
design meets the performance requirements, then a prototype of the selected medical device 
design may be built and tested. In addition, the visual result on the screen display 190 can be 

20 used by a physician to aid in the selection of various versions (e.g., sizes) of a given medical 
device design. For example, prior to a procedure for placement of a TPEG in a patient's aorta, 
the physician may first virtually test the performance of various TPEG designs or various 
versions of a single TPEG design prior to the procedure. To accomplish this, the physician 
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would obtain volumetric data from the patient's aorta by any of the various methods discussed 
above and input that data into an embodiment of a system 105 (in Figure 1) for analyzing the use 
of a medical device. The same or similar type of volumetric and materials data for a version of 
TPEG design to be tested is also loaded into the system 105. Note that it may be possible to load 
5 volumetric data from several anatomical features and versions of TPEG designs to be analyzed at 
one time, and then for the physician to choose which two to test together at a later time. Once 
the input data is loaded into the system 105, the visual result of the analysis of the 
Stress/Strain/Deformation Analyzer 160 is viewed by the physician on the screen display 190 
and based on those results, the physician determines whether the TPEG version tested meets, 

10 exceeds, or falls short of the clinical requirements of the patient. 

If the version of the TPEG which was virtually tested by the system 105 falls short of the 
clinical requirements of the patient, another version may be tested and so on until an appropriate 
design is identified. The physician may then begin the actual procedure on the patient with the 
appropriate TPEG design version. The system 105 may be configured to display the 

1 5 performance of a given TPEG design version with regard to long term structural integrity, 

prevention of perigraft leaks or sealing function, the general sizing of the TPEG with respect to 
the patient's aorta and the like. With regard to testing of the long term durability or structural 
integrity of the TPEG or other medical device design, the system 105 has great utility. 
Specifically, system 105 has the ability, assuming the use of sufficiently powerful CPUs, to 

20 recreate large numbers of cyclic expansions and contractions in a short period of time. For 
example, as discussed above, the vascular system of a patient is constantly expanding and 
contracting as a result of dynamic or static pressure gradients within the vasculature from the 
patient's beating heart. These expansions and contractions can put stresses, strains, and 
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deformations on intracorporeal medical devices, such as TPEG, which over time can lead to 
failure of the device. System 105 would give the physician the ability to quickly test a chosen 
TPEG design in a virtual model of the patient's expanding and contracting aorta for an amount of 
cycles that would equal or exceed the amount of cycles that would be expected in the patient's 
5 lifetime to determine the long term safety and efficacy of the design choice. Of course, a similar 
time compressed analysis could be used for any other type of medical device in any other part of 
a patient's body. Accordingly, if the invention is used as a preprocedure tool, physicians may 
analyze the use of various TPEG embodiments and select those that meet their performance 
requirements thereby allowing the physicians to select the best medical devices, such as the best 

1 0 TPEGs for treating their patients with aneurysm. 

Because of the computing resources needed by FEA software tools, they are generally run 
on Silicon Graphics or other UNIX computer systems. The Mesh Generator, 
Stress/Strain/Deformation Analyzer, and the visualization of the stresses, strains, and 
deformations on the candidate TPEG have been run on a Silicon Graphics (R12000) machine 

1 5 with 640MB of memory. 

Modifications to DYNA3D or NIKE3D 

In one embodiment, NIKE3D and DYN/A3D were used and modified to implement the 
features of the present invention (TPEG design was analyzed). In determining the required 
material model, an exemplary material model (herein called TPEG material model (W)) was used 

20 to accommodate a strain energy density of the form: 



19 



24641-1070 



W = a (7-3) + a (7 -3) + a (7-3) +a (7-3)(7-3) + a (1-3) + 

10 V 1 J 01 V 2 J 20 V 1 ' 11 V 1 JK 2 J 02 v 2 J 

a (7-3) + a (7 -3) \l -3) + a (7 -3)(7 -3) 2 + a (7-3) 3 + 

30 V 1 ' 21 v 1 ' v 2 7 12 v 1 /v 2 7 03 v 2 ' 

with 7T=2(a +a )/(l-2v) 
v io or v J 

5 

where 

a., are material parameters; 
v is Poisson's ratio; 

K is the bulk modulus given as a function of Poisson's ratio; and 
10 7 , 7 2 , and 1^ are the first, second, and third invariants of the right Cauchy-Green 

strain tensor, respectively. 

The TPEG material model (W), discussed above, was derived from a doctoral thesis, 
which discusses the stress in abdominal aortic aneurysm. (See Madhavan Lakshmiraghavan, 
Mechanical Wall Stress in Abdominal Aortic Aneurysm: Towards Development of a Clinical 

1 5 Tool to Predict Aneurysm Rupture (1 998) (unpublished Ph.D. dissertation, University of 
Pittsburgh which is hereby incorporated herein in its entirety). 

Other articles discussing a hyperelastic material, linear elastic, and non-linear elastic 
models of the aortic walls may also be used to derive a material model as exemplified above and 
other applications of the virtual prototyping system 105 (in Figure 1). (See M. L. Raghavan et 

20 al., Ex Vivo Biomechanical Behavior of Abdominal Aortic Aneurysm: Assessment Using a New 
Mathematical Model, 24 Annals of Biomedical Engineering 573-582 (1996); David A. Vorp. Et 
ah, Finite Element Analysis of the Effect of Diameter and Asymmetry on the Wall Stress 
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Distribution in Abdominal Aortic Aneurysm, 35 BED (Bioengineering Conference ASME 1 997) 
33-34 (1997), both of which are incorporated by reference herein in their entirety). 

Modifications to NIKE3D 

5 NIKE3D has an existing material model, number 15, which is a three-dimensional 

continuum hyperelastic material that uses a strain energy density function of the form: 

W = A(I -3) + B(7 -3) + ViKQaff) 2 



10 



with £"= 4(A+BXl+v) 
(3-Gv) 



where 

A and B are material parameters; 
v is Poisson's ratio; 

Kis the bulk modulus given as a function of Poisson's ratio; 
15 J and I are the first and second invariants of the right Cauchy-Green strain 

tensor, respectively; and 

9 is the current volume of the element divided by the undeformed volume. 
Using the material model 15 as the framework, the material model 15 is modified to 
implement the TPEG Material Model "W" listed above. This entails ensuring that variables are 
20 accordingly updated or modified in the source code to capture the information required by the 
TPEG Material Model. Material model 15 was chosen from the NIKE3D models because it 
involves the least amount of code modification to implement the features of the present 
invention. 
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Implementation of the TPEG Material Model in NIKE3D 

To implement the features in accordance with the present invention, two N1KE3D 
subroutines, wevahf and printm.f, were modified. 

The following modifications were made to NIKE3D subroutine weval.f: 
5 a) Ten material parameters (a 10? a^ a 11? a^ a 30 a 12 , a 03 K) were read instead of 

three (A, B, mdK). 

b) The calculation of K was changed from K = 4(A+B)(l+v)) / (3 - 6v) to 
^=2(a 10 + a 01 )/(l~2v) 

c) The calculation of dW was changed from 5W = A to 
10 dl, dl x 

SW = a 10 + 2a 20 (/ r 3) + a n (/ 2 -3) + 3a 30 (/ r 3) 2 + 2a 21 (/ r 3)(/ 2 -3) + a l2 (/ 2 -3) 2 
dl x 

15 d) The calculation of dW was changed from dW = B to 

8I 2 8I 2 

dW_ = ao, + a n (7 r 3) + 23^-3) + ^(I^f + 2a !2 (7 r 3)(/ 2 -3) + 3ao 3 (7 2 -3) 2 



20 



a/ 2 



e) The higher derivatives of W with respect to 7, and 7 2 were changed from zero to 
= 20,0 + 6830(7,-3) + 232,(7,-3), 

57, 2 

25 3*W = 2ao 2 + 2a, 2 (7,-3) + 6803(7,-3), and 

dL 2 



30 



= a„ + 232,(7,-3) + 2a 12 (7 2 -3) 

dl x dl 2 
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f) The derivatives with respect to I 3 were changed from 5W = K (\nl 3 1 1 3 ) to 



81, 



gW = K(I 3 - 1) and from dV = K ((1 - ln/ 3 )//j 2 ) to <fW= K 
31 dl 2 dlf 

3 3 J 

g) When a completely incompressible material (7=1) is specified by setting the 
10 augmented Lagrangian flag to true, the derivatives with respect to are left in the 

log form. The log form shows substantially faster convergence and better stability 
for completely incompressible materials. 
The NIKE3D subroutine printm.f was modified to print out all nine a material 

parameters to the material description in the high-speed printout file. 
1 5 Invocation of the Modified NIKE3D TPEG Material Model 

The TPEG material model (W) (i.e., the modified NIKE3D Material Model 15) is 
invoked in NIKE3D using the input data format shown in Table I. Poisson's ratio is kept as the 
third parameter to maintain compatibility with models using the original NIKE3D hyperelastic 
model. The documentation for NIKE3D, and the TRUEGRID Mesh Generator, provides an 
20 input format list for Material Model 1 5 similar to Table I given below, with A, B, and v all 
defined on card 3 (it should be understood that the "card" represents lines of input data). The 
original NIKE3D code, however, reads A from columns 1-10 card 3, B from columns 1-10 of 
card 4, and v from columns 1-10 of card 5. This format was changed to comply with the 
NIKE3D manual and the format in Table I in the modified weval.f and printm.f subroutines. 

25 
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Table I - 

Input parameters format for the modified NIKE3D material model (TPEG material model) 


Card 


Columns 


Description 


Format 


1 


1-5 


Material ID number 


15 


1 


6-10 


Material type (use 15) 


15 


1 


11-20 


Density 


E 10.0 


1 


21-25 


Element class (not used) 


15 


1 


26-35 


Reference temperature (not used) 


E 10.0 


1 


36-45 


Rayleigh damping parameter alpha 


E 10.0 


1 


46-55 


Rayleigh damping parameter beta 


E 10.0 


2 


1-72 


Material title 


12A6 


3 


1-10 


a 

10 


E 10.0 


3 


11-20 


a 

01 


E 10.0 


3 


21-30 


Poisson's ratio 


E 10.0 


3 


31-40 


a 

20 


T7 1 A A 

Jb 10.0 


3 


A 1 CA 

41-50 


a 

ii 


U 1 A A 

E 1U.0 


D 


J 1-00 


a 

02 


"p i n a 


3 


61-70 


a 

30 


E 10.0 


3 


71-80 


a 

21 


E 10.0 


4 


M0 


a 

12 


E 10.0 


4 


11-20 


a 

03 


E 10.0 
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5-7 


All 


Blank 




8 


1-10 


Augmented Lagrangian flag 

.EQ.l : active, enforce compressibility with augmented 
Lagrangian iteration 


E10.0 


8 


11-20 


Convergence tolerance for augmented Lagrangian 
iteration 

.GT.0.0: converged when volume strain norm < TOL 
(tolerance) 

.LT.0.0: augment exactly - TOL times 


E 10.0 



The format column specifies the expected data type. For example, a format of "I" means 
that an integer is expected ("15" means integer with 5 positions), "E" means a real numeric value, 
and "A" means character data type. 

Modifications to DYNA3D 

DYNA3D has an existing material model number 27, which is a three-dimensional 
continuum hyperelastic material that uses a strain energy density function of the form 

W - A(/-3) + B(J -3) + C(/ 2 - 3) + D(J - 3) 2 

with C - V% A + B 
and 

D = A(5v-2) + B(llv-5) 
2-4v 

where: 
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A and B are material parameters; 
v is Poisson's ratio; and 

/ ,1^, and J are the first, second, and third invariants of the right Cauchy-Green 

strain tensor, respectively. 
5 The material model 27 may be modified to implement the TPEG Material Model (W)). 

This also entails ensuring that variables are accordingly updated or modified in the source code 
to capture the information for the TPEG material model (W). 

Implementation of the TPEG Material Model in DYNA3D 

10 To implement the features in accordance with the present invention, two DYNA3D 

-2 

subroutines, f3dm27.f and printm.f, were modified. The "C(I - 1)" term was left in the 

modified material model since without it, the explicit time integrator becomes unstable very 
easily. This term only significantly changes the result when the material undergoes significant 
change in volume. If v « 0.5, the material behaves in a nearly incompressible matter, in this case 
15 D is much larger than C 7 and the inclusion of C has little to no effect on the final result. 

The following modifications were made to DYNA3D subroutine £3dm27.f: 

a) Ten material parameters (a 10 , a^ a^ a 11? a^ ^ a^, a 12 , a^ K) were read instead of four 
(A, B, C, andD). 

20 b) The calculation of D was changed from D = (A(5v - 2) + B(l 1 v - 5)) / (2 - 4v) to 

D = (a +a )/(l-2v) 

c) The computation for I ? and I were added. 
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d) The calculation of dW was changed from 6W = A to 
dl dl 



1 

3W =a +2a (7-3) + a (7 -3) + 3a (7-3) 2 + 2a (7-3)(7-3) + a (7-3) 2 . 

10 20 v 1 7 II x 2 7 30 v 1 7 21 v 1 7V 2 7 12 v 2 7 

57 

1 



e) The calculation of dW was changed from dW = B to 



5/ dl 

2 2 

10 aw =a +a (7-3) + 2a (7-3) + a (7-3) 2 + 2a (/ -3)(7 -3) + 3a (7-3)1 

01 11 V 1 7 02 x 2 7 21 v 1 7 12 v 1 7V 2 7 03 v 2 7 

dl 

2 



_3 

f) The calculation of 3W = 2D(7 -1) - 2C(7 - 1) remains unchanged, however, 
15 81 

3 

the value of D has changed. 
The DYNA3D subroutine printm.f was modified to correctly output the hyperelastic 
20 material constants to the resulting high-speed printout file. 



Invocation of the Modified DYNA3D Material Model (TPEG Material Model) 

The TPEG material model (i.e., the modified DYNA3D material model 27) is invoked in 
DYNA3D using the input data format shown in Table II. Poisson's ratio is kept as the third 
25 parameter to maintain compatibility with models using the original DYNA3D hyperelastic 
model. 
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Table II - 

Input parameters format for the modified DYNA3D material model (TPEG material model) 


Card 


Columns 


Description 


Format 


1 


1-5 


Material ID number 


15 


1 


6-10 


Material type (use 15) 


15 


1 


11-20 


Density 


E 10.0 




21-25 


Element class (not used) 


15 


1 


26-35 


Reference temperature (not used) 


E 10.0 


1 


36-45 


Rayleigh damping parameter alpha 


E 10.0 


1 


46-55 


Rayleigh damping parameter beta 


E 10.0 


2 


1-72 


Material title 


12A6 


3 


1-10 


a 

10 


E 10.0 


3 


11-20 


a 

01 


E 10.0 


3 


21-30 


Poisson's ratio 


E 10.0 


J 


■}1 ACS 


a 

20 


Finn 


"J 


1 - JU 


a 

ii 


Finn 


3 


51-60 

± \J\J 


o 
<x 

02 


E 10.0 


3 


61-70 


a 

30 


E10.0 


3 


71-80 


a 

21 


E 10.0 


4 


1-10 


a 

12 


E10.0 


4 


11-20 


a 

03 


E 10.0 
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5-7 


All 


Blank 





Reading the doctoral thesis mentioned above, the appropriate values of input parameters 
may accordingly be provided as input to the Stress/Strain/Deformation Analyzer {see Madhavan 
Lakshmiraghavan, Mechanical Wall Stress in Abdominal Aortic Aneurysm: Towards 
5 Development of a Clinical Tool to Predict Aneurysm Rupture (1998) (unpublished Ph.D. 
dissertation, University of Pittsburgh). 

TRUE GRID Command File 

Figures 5A through 5M contain a command file that is an exemplary file read by 
10 TRUEGRID to implement the features of the present invention (e.g., for stent design). This 

exemplary command file illustrates a component-level analysis of a stent, without the option for 
simulating deployment into CT-based anatomy (isim mode=6, not present in the command file). 

TRUEGRID, in its basic form, is not only a Mesh Generator, but is also a format 
generator. It outputs data in a certain format, which are then read by NIKE3D and/or DYNA3D. 
1 5 The invention utilizes both TRUEGROD's capability as a Mesh Generator and an output 
generator to create an output file (e.g., Tables I and II discussed above), containing the 
appropriate values that would be read by NIKE3D and DYNA3D, respectively. The outputs 
created by TRUEGRID may be created by other means, e.g., by other Mesh Generator software 
or proprietary software. 
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The command file (contained in Figures 5 A-5M) contains the parameters and the 
instructions that are read by TRUEGRID to generate the mesh and the output file(s), which are 
read by DYNA3D and/or NIKE3D. 

The line numbers at the start of each line are only added to facilitate reference to 
5 particular lines in the command file and are not part of the command file. Text after the "c " are 
ignored by TRUEGRID (comments). To take advantage of the capabilities of TRUEGRID, the 
command file contains various parameters that help developers customize their simulation and/or 
Stress/Strain/Deformation analysis. Mesh generating tools, such as TRUEGRID, in the non- 
interactive mode, generally require that command files or similar files be created to enable them 

10 to generate finite element models. In the interactive mode, a finite element model may be 

created by a medical device designer (e.g., TPEG designer) using the options available in the 
interactive mode of TRUEGRID. 

Referring to Figure 5 A, the inike parameter (lines 5 and 21) tells TRUEGRID that the 
output file is to be read by a NIKE3D Stress/Strain/Deformation Analyzer. The command file 

1 5 also tells TRUEGRID that the stent to be modeled is a full 3-segment stent design (line 6 and 
22), the model is a full 360 degree model of a stent (lines 6 and 23), to model the stress on the 
initial expansion of the stent in vivo (lines 16 and 24), and to refine the elements by 2 in each 
direction of the cross section (lines 1 8 and 25). (Crowns can be a pointed or barbed portion of a 
stent - see lines 7 through 9). The command file thus enables TRUEGRID to generate a mesh 

20 and a model of a stent subjected to various component-level in vitro tests such as radial force and 
predelivery compression. Simulation of these tests enables a designer to refine and optimize the 
stent design for its intended application (e.g. as component of a TPEG or for treating occlusive 
disease). 
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TRUEGRID can also act like an interpreter. It reads the information contained in the 
command file, and interprets and processes the lines accordingly. For example, the text after the 
word "para" or "parameter" are parameters read by TRUEGRID. These terms indicate the value 
or the formula that should be used by TRUEGRID. For example, line 21 denotes that the 
5 parameter inike contains the initial value 1 . 

Line 46 in Figure 5B means that the value of the parameter dCIA3 contains the value 0.0. 
Line 138 in Figure 5D indicates that the initial value of the parameter rocompcyl is the 
value evaluated by the formula "[0.95*(min(%RCyl3 ? %RCyl6,%RCyll2_l ? %RCyll2_2)- 
%RW6)." TRUEGRID understands that the min function has to be evaluated. The min function 
10 compares the value contained in each variable, in this case, contained in RCyl3 (e.g., contains 
1), RCyl6 (contains 0.005), RCyll2_l (contains 0.987), and RCyll2_2(contains 0.0002), and 
returns the content of the variable, which holds the least value — 0.0002 (value contained in 
Rcyll2_2). Assuming the variable RW6 contains the value 0.18, TRUEGRID then evaluates the 
rocompcly variable to contain 0.95 * 0.0002 - 0.18, which equals to negative 0.17981. This 
1 5 value is thus the initial value of rocompcyl when initially processed and read by TRUEGRID. 

Embodiments of the invention can simulate various phases of TPEG use. For example, it 
calculates the stresses, strains, and deformations on the TPEG when it is compressed then 
decompressed for deployment, when the TPEG is compressed into the catheter for deployment, 
when the TPEG expands, and the like. 
20 Referring to line 432, in Figure 5L, the term "include" indicates to TRUEGRID that 

when the condition as defined in line 43 1 is met, the istent.mts nike solid file is read. The 
contents of this include file could be added in the command file itself. For flexibility and 
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readability, however, they were placed in a separate file. Programmers typically use include 
files, such as done in C or C++, for code control and ease of maintenance 

Figures 6A-6F contain an exemplary text of a command file called "seal.run" (line 2) read 
by TRUEGRID for a simulated TPEG graft deployment in a proximal aortic neck to generate a 
5 mesh incorporating both an anatomical feature and medical device and to output files that are 
read by a Stress/Strain/Deformation Analyzer. 

Figures 7A-7C is an exemplary include file, called "tpeg.part_ct_aorta3," used by 
"seal.run" command file listed in Figures 6A-6F. See line 217 of Figure 6F. This file contains 
the commands which read in surfaces created by the Geometry Generator 120 from CT data for 
1 0 the aorta and builds the mesh for the vessel. 

Figures 8 A-8L is another exemplary command file read by TRUEGRID used in the 
virtual prototyping system of the present invention for simulating stent deployment into an 
anatomy from CT-data, as opposed to a stent graft. The stent could be a part of a stent graft, 
could be intended for use to treat occlusive disease in the vasculature, or could even be used for 
1 5 nonvascular application, such as an esophageal stent. 

The files listed in Figures 5A-5M, 6A-6F, 7A-7C, and 8A-8L are written to be read by 
TRUEGRID. Variations on such files are expected depending on the Mesh Generator 130 
deployed in the system. 

Figure 9A illustrates a flow chart, which sets forth the basic components of an 
20 embodiment of the inventive system and process in accordance with the present invention. In 
particular, this figure illustrates how to develop better-designed TPEGs. The steps illustrated 
may of course be utilized for developing other medical devices, other than TPEGs. 
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To start, a TPEG designer first determines, in box 905 A, the performance requirements 
desired, such as to secure an optimal structural integrity of the TPEG, to avoid potential health 
risks such as ruptures and endoleaks, or to have a smaller TPEG packaging. 3D volumetric data 
of the anatomy desired, for example, in this case a blood vessel, is then acquired at box 91 OA, 
5 using CT or MRI scanners. Alternatively, if 3D volumetric data are already available, such 
acquisition may be skipped and such 3D volumetric data be obtained from the archive. 

It should be noted here that the "anatomy" desired, which defines the embodiment in 
which a medical device is to be tested, is not necessarily limited to a patient's body. For 
example, embodiments of the present invention could be used to obtain test results for medical 

10 device performance in a wide variety of in vitro tests, some of which may be necessary or 

desirable for Food and Drug Administration (FDA) approval of the medical device in question. 
Various forms of in vitro failure mode testing such on tensile pull testing and the like could be 
performed by an embodiment of the invention and allow the tester to easily vary test parameters, 
device design, and test frequency to quickly obtain the desired test results. In addition, 

1 5 volumetric anatomical data for animals could be used to simulate animal testing that is necessary 
or desirable for FDA approval of a medical device. This may be of particular importance for a 
medical device design, which seeks to establish equivalence with an existing approved product 
which has been previously tested in animal studies. 

The geometry generator (120 in Fig. 1) then generates a blood vessel geometric model in 

20 box 920A. As discussed above, the blood vessel geometric model may be an actual idealized or 
in vitro model. If the geometry generator is an embodiment where surface points are first 
extracted, a CAD system may then be used to generate such geometric model. 
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Next, a candidate TPEG model or design, which is obtained typically from a model 
created using a CAD software, is selected or modeled by the TPEG designer (step 925 A). The 
Mesh Generator (130 in Fig. 1) then generates a mesh model incorporating both the blood vessel 
and the TPEG (93 OA). A TPEG designer then determines the material properties of the candidate 
5 TPEG model and the blood vessel at step 935A. The material properties may also have been 
assigned by the TPEG designer during the previous step (i.e., the generation of the mesh model). 
Using a Stress/Strain/Deformation Analyzer (160 in Fig. 1), assuming that the load (150 in Fig. 
1) and the Materials Model (170 in Fig. 1) are available to the Stress/Strain/ Deformation 
Analyzer for input, a TPEG designer then simulates the candidate TPEG design behavior in a 

10 stress/strain/deformation analysis (at step 940A) to determine if the candidate TPEG meets the 
performance requirements. 

If the candidate TPEG does not meet the performance requirements, a "no" outcome at 
decision box 955A, the TPEG designer chooses another TPEG design or model at step 980A, and 
repeats the steps as shown by the arrow to box 925 A. If it, however, meets the target 

15 performance requirements, a "yes" outcome at decision box 955 A, a prototype is then fabricated 
based on the candidate TPEG model and design at step 960A. The fabricated prototype is then 
subjected to testing, e.g., animal testing or clinical testing, at step 965 A. If the fabricated 
prototype meets the target performance requirements, the candidate TPEG model thus is a final 
design and may be used to produce other TPEGs. 

20 If the fabricated prototype, however, does not meet the performance requirements, a "no" 

outcome at decision box 970A, the TPEG designer modifies the TPEG design or selects a new 
TPEG design, and repeats the steps as shown with the arrow to box 925 A. If necessary, the 
process is repeated several times until the performance requirements and the final design is 
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obtained. A benefit of the invention is to reduce the number of "no" outcome at decision box 
970A compared to a development process which uses only hardware prototypes for design 
verification. 

As discussed above, a proposed TPEG model may be evaluated against a number of 
5 anatomical features to determine the suitable range of conditions of an applicable TPEG model 
(e.g., size). Similarly, a set of anatomical features may be evaluated against a number of TPEG 
models to determine the type of suitable TPEG model for such set of anatomical feature. 
Furthermore, an analysis of the stresses, strains, and deformations may be conducted on the 
medical device without interaction to certain anatomical features. 

10 Figure 9B, is similar to figure 9A except for the additional step (box 942B) of displaying 

the visual simulation of the stresses and strains on the TPEG. The display of the simulation is 
typically employed using the Visualization tool (180 in Figure 1), which in the preferred 
embodiment is the GRIZ software. 

Visual display of the simulation is not necessary because a reading of the numerical 

1 5 representation of the stresses, strains, and deformation on the TPEG may guide a TPEG designer 
whether the performance requirements are met. However, visual display is often desirable 
because a visual representation of the stresses and strains, for example, red hot spots on the 
visual TPEG model can be easier to understand than mere numerical representations. 

Figure 10 is similar to Figure 9 A and illustrates a process to develop better-designed 

20 medical devices using in vitro features. In the first step as shown in 1005, a medical device 
designer, determines the performance requirements. The next step is to generate a geometry 
model of the in vitro model, step 102 OA, (e.g., latex tube to represent an artery), using software 
tools, such as a CAD software or even TRUEGRID. The steps are then similar to those 
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illustrated in Figure 9A, In another embodiment, the in vitro model such as a latex tube may be 
scanned to obtain a 3D volumetric data. Such acquired 3D volumetric data may also be modified 
by the medical device designer. 

In another embodiment not shown, only the medical device model is analyzed absent the 
5 anatomical feature or in vitro model. The operations shown in Figure 10 would be implemented, 
without the operation of generating blood vessel geometric model (step 1020A) and the analysis 
would only be performed on the geometric model of the candidate medical device or a potion of 
it. Material properties and load information pertinent only to the medical device are generally 
used in the analysis process. 

10 Figure 1 1 contains steps similar to those illustrated in Figure 9A. Figure 1 1 illustrates an 

embodiment of the present invention as a preprocedure planning tool, for example, to guide a 
physician in deciding which particular TPEG to implant in a patient. 

To start, a physician first determines, in box 1 1 05, the surgical or interventional 
procedure objectives, typically, to ensure robust sealing and structural integrity of the TPEG in 

1 5 vivo for a particular patient. The physician then obtains 3D volumetric data of the potential site 
of the TPEG, e.g., the abdominal aorta, at step 1 1 10. The Geometry Generator (120 in Fig. 1) 
then extracts the surface points from the 3D volumetric data acquired in step 1115. Based on the 
surface points extracted, a blood vessel geometric model is created 1 120. 

Next, a candidate TPEG, which is obtained typically from a model created using a CAD 

20 software, is selected by the physician (step 1 125). (TPEG models may be created in advance and 
stored in a library in the system. At this point, the physician is determining which available 
TPEG design is best suited for that patient or individual). The Mesh Generator (130 in Fig. 1) 
then generates a mesh model incorporating both the blood vessel and the selected TPEG. A 
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physician may then identify the material properties of the candidate TPEG and the blood vessel 
at step 1 135. The material properties may have also been assigned during the previous step (i.e., 
the generation of the mesh model). Using a Stress/Strain/Deformation Analyzer (160 in Fig. 1), 
assuming that the load (150 in Fig. 1) and the materials model (170 in Fig. 1) are available to the 
5 Stress/Strain/Deformation Analyzer for input, a physician may then run the candidate TPEG to a 
stress/strain/deformation analysis (at step 640) to determine if the candidate TPEG meets the 
surgical objectives. 

If the candidate TPEG does not meet the procedural objectives, a "no" outcome at 
decision box 1 155, a physician may decide to change the TPEG to be used in the procedure at 

1 0 step 1 1 80 and repeat the process as shown by the arrow to box 1 125. Based on the physician's 
judgment, if the candidate TPEG does meet the prcedural objectives, a "yes" outcome at decision 
box 655, the physician then may decide whether to proceed with the planned TPEG implant 
procedure or not, at step 1 160. 

Figure 12 contains a representation of one simulation display of a cutaway lateral 

1 5 view of a vascular stent in the infrarenal aorta just proximal to an abdominal aneurysm. 
Using the system as described above, several displays may be presented to the user 
showing the progressive stent expansion and contact with the luminal surface of the 
vessel. The system may be also be used such that the visualization module displays the 
medical device and the anatomical feature in color, with colors and their gradients 

20 representing the various stresses, strains, and deformations on the medical device and the 
anatomical feature. Other views, such as a proximal view, may also be used in 
simulation. Figure 13 is a block diagram of an exemplary computer 1300 such as might 



37 



24641-1070 




comprise any of the computers containing a Geometry Generator 120, a Mesh Generator 130, a 
Stress/Strain/Deformation Analyzer 160, and a Visualization tool 180. Each computer 1300 
operates under control of a central processor unit (CPU) 1302, such as a high-end microprocessor, 
e.g., typically found in Silicon Graphics workstation, and associated integrated circuit chips, A 
5 computer user can input commands and data from a keyboard and mouse 1312 and can view inputs 
and computer output at a display 1310. The display is typically a video monitor or flat panel 
display device. The computer 1300 also includes a direct access storage device (DASD) 1304, such 
as a fixed hard disk drive. The memory 1306 typically comprises volatile semiconductor random 
access memory (RAM). Each computer preferably includes a program product reader 1314 that 

1 0 accepts a program product storage device 1316, from which the program product reader can read 
data (and to which it can optionally write data). The program product reader can comprise, for 
example, a disk drive, and the program product storage device can comprise removable storage 
media such as a floppy disk, an optical CD-ROM disc, a CD-R disc, a CD-RW disc, DVD disk, or 
the like. In the preferred embodiment, each computer 1300 can communicate with the other 

1 5 connected computers over the network 1 320 through a network interface 1 308 that enables 

communication over a connection 1318 between the network and the computer. This facilitates 
having each separate system as illustrated in Figure 1, provide inputs and outputs to the other 
components in the system. 

The CPU 1302 operates under control of programming steps that are temporarily stored in 

20 the memory 1 306 of the computer 1 300. When the programming steps are executed, the pertinent 
system component performs its functions. Thus, the programming steps implement the 
functionality of the system components illustrated in the figures. The programming steps can be 
received from the DASD 1304, through the program product 1316, or through the network 
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connection 1318. The storage drive 1304 can receive a program product ? read programming steps 
recorded thereon, and transfer the programming steps into the memory 1306 for execution by the 
CPU 1302. As noted above, the program product storage device can comprise any one of multiple 
removable media having recorded computer-readable instructions, including magnetic floppy disks, 
5 CD-ROM, and DVD storage discs. Other suitable program product storage devices can include 
magnetic tape and semiconductor memory chips. In this way, the processing steps necessary for 
operation in accordance with the invention can be embodied on a program product. 

Alternatively, the program steps can be received into the operating memory 1306 over the 
network 1318. In the network method, the computer receives data including program steps into the 

10 memory 1306 through the network interface 1308 after network communication has been 

established over the network connection 1318. The program steps are then executed by the CPU 
1302 to implement the processing of the present invention. 

Although the present invention is implemented on UNIX workstations, typical personal 
computers could likely be adopted to perform these functions in the future. 

15 It should be understood that all of the computers of the systems embodying the various 

systems illustrated in Figure 1 , preferably have a construction similar to that shown in Figure 1 3, so 
that details described with respect to the Figure 13 computer 1300 will be understood to apply to all 
computers or components of the system. Any of the computers can have an alternative 
construction, so long as they have sufficient resources and processing power to handle finite 

20 element analyses and other functions in accordance with the present invention. 

Those skilled in the art will recognize that variations in the steps, as well as the order of 
execution, may be done and still make the various embodiments of the invention operate. 
Furthermore, one skilled in the art will realize that although the examples described herein 
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generally refer to TPEGs, other medical devices may be designed in accordance with the present 
invention. 

In addition, although the modules of the system 105 (Figure 1), the Geometry Generator, 
the Mesh Generator, Stress/Strain/Deformation Analyzer, and the Visualization module, are 
5 shown in different boxes, depending on the software tools utilized their functions may with each 
other. Some functions, for example, that are done by one module, e.g., the Mesh Generator, 
TRUEGRID, thus, may also be done by the Geometry Generator, MIMICS, or vice versa. 

Embodiments of the present invention have been described above so that an 
understanding of the present invention can be conveyed. There are, however, many alternative 
10 software programs available or able to be written that would embody the functions of the present 
invention, and thus, may be used accordingly. The present invention should therefore not be 
seen as limited to the particular embodiments described herein, but rather, it should be 
understood that the present invention has wide applicability with respect to medical device 
design generally. All modifications, variations, or equivalent arrangements and implementations 
1 5 that are within the scope of the attached claims should therefore be considered within the scope 
of the invention. 
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is claimed is: 

A system for analyzing the use of medical devices comprising: 

a) geometry generator that receives three-dimensional volumetric data of at least one 
anatomical feature and generates a geometric model of said anatomical feature; 

b) mesh generator that receives the said geometric model of said anatomical feature and 
the geometric model of a medical device, and generates a finite element model or mesh 
incorporating both said anatomical feature and said medical device; and 

c) stress/strain/deformation analyzer that receives said mesh incorporating both said 
anatomical feature and said medical device, materials properties of said anatomical 
feature and said medical device, and load on said anatomical feature and/orsaid medical 
device, and simulates stresses, strains, and deformations of said medical device. 

A system as defined in claim 1 where said geometric model of said anatomical feature is 
an idealized geometric model. 

A system as defined in claim 1 where said three-dimensional volumetric data are acquired 
via CT scan. 

A system as defined in claim 1 where said three-dimensional volumetric data are acquired 
via MRL 

A system as defined in claim 1 where said geometric model of a said medical device is 
for an endovascular prosthesis. 

A system as defined in claim 5 where said endovascular prosthesis is a transluminally 
placed endovascular graft. 
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7. A system as defined in claim 5 where said endovascular prosthesis is a cardiovascular 
stent device. 

8. A system as defined in claim 1 where said geometry generator is MIMICS. 

9. A system as defined in claim 1 where said mesh generator is TRUEGRID. 

10. A system as defined in claim 1 where said stress/strain/deformation analyzer is 
DYNA3D. 

11. A system as defined in claim 1 where said stress/strain/deformation analyzer is NIKE3D. 

12. A system as defined in claim 10 where said DYNA3D is modified to accommodate a 
strain energy density of the form: 

W =a (7-3) + a (7 -3) + a (7-3) 2 + a (7-3)(/-3) + a (I-3) 2 + 

10 v 1 7 01 v 2 7 20 v 1 7 ll v 1 7V 2 7 02 v 2 7 

a (7 -3) + a (7-3) 2 (7-3) + a (7 -3)(7 -3) 2 + a (7-3) 3 + 

30 v 1 7 21 v 1 7 v 2 7 12 v 1 7V 2 7 03 v 2 7 

with i:=2(a +a )/(l-2v) 

v jo or v 7 

where 

a., are material parameters; 
v is Poisson's ratio; 

K is the bulk modulus given as a function of Poisson's ratio; and 

7 , 7 2? and 1^ are the first, second, and third invariants of the right Cauchy-Green 

strain tensor, respectively. 
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13. A system as defined in claim 1 1 where said NIKE3D is modified to accommodate a strain 
energy density of the form: 

W = a (7-3) + a (7-3) + a (7-3) 2 + a (7-3)(7-3) + a (I-3) 2 + 

10 v 1 7 01 v 2 7 20 v 1 7 \\ x 1 7V 2 7 02 v 2 7 

a (7-3) + a (7 -3) 2 (7 -3) + a (7 -3)(7 -3) 2 + a (7-3) 3 + 

30 V 1 7 21 V 1 7 V 2 7 12 V 1 7V 2 7 03 v 2 7 

5 YiK(I-Vf 

with 7T = 2(a +a )/(l-2v) 

V 10 01 7 v 7 

where 

a.^ are material parameters; 

10 v is Poisson's ratio; 

K is the bulk modulus given as a function of Poisson's ratio; and 

7 , 7 2 , and 1^ are the first, second, and third invariants of the right Cauchy-Green 

strain tensor, respectively. 

14. A system as defined in claim 1 further comprising visualization tool that receives said 
1 5 stresses and strains on said medical device and anatomical feature and displays said 

stresses and strains on said medical device via visual representation. 

15. A system as defined in claim 14 where said visualization tool is GRIZ. 

16. A system for analyzing the use of a medical device comprising: 

a) geometry generator that receives three-dimensional volumetric data of at least one 
20 anatomical feature of a particular individual and generates a geometric model of said 

anatomical feature; 
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b) mesh generator that receives the said geometric model of said anatomical feature and 
the geometric model of a medical device, and generates a finite element model or mesh 
incorporating both said anatomical feature and said medical device; and 

c) stress/strain/deformation analyzer that receives said mesh incorporating both said 
5 anatomical feature and said medical device, materials properties of said anatomical 

feature and said medical device, and load on said anatomical feature and/or said medical 
device, and simulates stresses, strains, and deformation of said medical device. 
17. A system as defined in claim 16 where said geometric model of said anatomical feature 
is an idealized geometric model. 
10 18. A system as defined in claim 16 where said three-dimensional volumetric data are 
acquired via CT scan. 

19. A system as defined in claim 16 where said three-dimensional volumetric data are 
acquired via MRL 

20. A system as defined in claim 16 where said geometric model of a said medical device is 
1 5 for an endovascular prosthesis. 

21 . A system as defined in claim 20 where said endovascular prosthesis is a transluminally 
placed endovascular graft. 

22. A system as defined in claim 20 where said endovascular prosthesis is a cardiovascular 
stent device. 

20 23. A system as defined in claim 16 where said geometry generator is MIMICS. 

24. A system as defined in claim 16 where said mesh generator is TRUEGRID. 

25. A system as defined in claim 16 where said stress/strain/deformation analyzer is 
DYNA3D. 
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26. A system as defined in claim 16 where said stress/strain/deformation analyzer is 
NIKE3D. 

27. A system as defined in claim 25 where said DYNA3D is modified to accommodate a 
strain energy density of the form: 

5 W =a (7-3) + a (7-3) + a (7-3) 2 + a (7-3)(7-3) + a (I -3) 2 + 

10 v \ 7 01 v 2 J 20 v 1 7 ll y 1 /v 2 J 02 v 2 J 

a (7-3) + a (7 -3) \l -3) + a (7 -3)(7 -3) 2 + a (7-3) 3 + 

30 v 1 ' 21 v 1 y v 2 y 12 v 1 yv 2 y 03 v 2 7 

with i^=2(a +a )/(l -2v) 
1 0 where 

a„ are material parameters; 
v is Poisson's ratio; 

X is the bulk modulus given as a function of Poisson's ratio; and 

7 , 7 2? and 7^ are the first, second, and third invariants of the right Cauchy-Green 

1 5 strain tensor, respectively. 
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28. A system as defined in claim 26 where said NIKE3D is modified to accommodate a strain 
energy density of the form: 

W =a (7-3) + a (7-3) + a (7-3) 2 + a (7-3)(7-3) + a (I-3) 2 + 

10 V 1 7 01 V 2 7 20 V 1 7 11 V 1 7V 2 7 02 V 2 7 

a (7 -3) + a (7 -3) 2 (7 -3) + a (7 -3)(7 -3) 2 + a (7-3) 3 + 

30 v 1 7 2\ K 1 7 v 2 7 12 x 1 7V 2 7 03 v 2 7 

5 ^(^-l) 2 

with 7T=2(a +a )/(l -2v) 

V 10 01 7 v 7 

where 

a are material parameters; 

10 v is Poisson's ratio; 

K is the bulk modulus given as a function of Poisson's ratio; and 

7 , 7 2? and 1^ are the first, second, and third invariants of the right Cauchy-Green 

strain tensor, respectively. 

29. A system as defined in claim 16 further comprising visualization tool that receives said 
1 5 stresses and strains on said medical device and anatomical feature and displays said 

stresses and strains on said medical device via visual representation. 

30. A system as defined in claim 29 where said visualization tool is GRIZ. 

31. A system for analyzing the use of medical device comprising: 

a) mesh generator that receives a geometric model of in vitro feature and a geometric 
20 model of a medical device, and generates a finite element model or mesh incorporating 

both said in vitro feature and said medical device; and 

b) stress/strain/defonnation analyzer that receives said mesh incorporating both said 
anatomical feature and said medical device, materials properties of said anatomical 
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feature and said medical device, and load on said anatomical feature and/or said medical 
device, and simulates stresses, strains, and deformations on said medical device. 
32. A system as defined in claim 3 1 where said in vitro feature is a geometric model of an 
idealized anatomical feature. 
5 33. A system as defined in claim 3 1 where said geometric model of said medical device is for 
an endovascular prosthesis. 

34. A system as defined in claim 33 where said endovascular prosthesis is a transluminally 
placed endovascular graft. 

35. A system as defined in claim 33 where said endovascular prosthesis is a cardiovascular 
1 0 stent device. 

36. A system as defined in claim 31 where said mesh generator is TRUEGRID. 

37. A system as defined in claim 31 where said stress/strain/deformation analyzer is 
DYNA3D. 

38. A system as defined in claim 31 where said stress/strain/deformation analyzer is 
15 NIKE3D. 
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39. A system as defined in claim 37 where said DYNA3D is modified to accommodate a 
strain energy density of the form: 

W =a (7 -3) + a (7 -3) + a (7-3) 2 + a (7-3)(7-3) + a (I-3) 2 + 

10 x I 7 01 v 2 7 20 v 1 7 ll v 1 /v 2 7 02 v 2 7 

a (7-3) + a (/ -3) \l -3) + a (7 -3)(7 -3) 2 + a (7-3) 3 + 

30 V 1 7 21 V 1 7 V 2 7 12 v 1 7V 2 7 03 V 2 7 

5 ^(/-l) 2 

with JC = 2(a +a )/(l -2v) 

v io or v 7 

where 

a are material parameters; 
10 v is Poisson' s ratio ; 

K is the bulk modulus given as a function of Poisson' s ratio; and 

7 i? 1 9 and 7 3 are the first, second, and third invariants of the right Cauchy-Green 

strain tensor, respectively. 

40. A system as defined in claim 38 where said NIKE3D is modified to accommodate a strain 
1 5 energy density of the form: 

W =a (7-3) + a (J -3) + a (7-3) 2 + a (7-3)(7-3) + a (I-3) 2 + 

10 v 1 7 01 v 2 7 20 v 1 7 ll v 1 7V 2 7 02 v 2 7 

a (7 -3) + a (7 -3) 2 (7 -3) + a (7 -3)(J -3) 2 + a (7-3)' + 

30 v 1 7 21 v 1 7 v 2 7 12 V 1 /x 2 7 03 v 2 7 

^(7-l) 2 

with *T = 2(a +a )/(l-2v) 

v 10 or v 7 



where 



a., are material parameters; 



v is Poisson' s ratio; 
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K is the bulk modulus given as a function of Poisson's ratio; and 

7 , 1 , and I are the first, second, and third invariants of the right Cauchy-Green 

strain tensor, respectively. 

41 . A system as defined in claim 31 further comprising visualization tool that receives said 
stresses and strains on said medical device and anatomical feature and displays said 
stresses and strains on said medical device via visual representation. 

42. A system as defined in claim 41 where said visualization tool is GRIZ. 

43. A system for analyzing the use a of medical device comprising: 

a) mesh generator that receives a geometric model of a medical device, and generates a 
finite element model or mesh of said medical device; and 

b) stress/strain/deformation nonlinear analyzer that receives said mesh, materials 
properties of said medical device, and load on said medical device, and simulates stresses, 
strains, and deformations on said medical device. 

44. A system as defined in claim 43 where said geometric model of said medical device is for 
an endo vascular prosthesis. 

45. A system as defined in claim 44 where said endo vascular prosthesis is a transluminally 
placed endovascular graft. 

46. A system as defined in claim 44 where said endovascular prosthesis is a cardiovascular 
stent device. 

47. A system as defined in claim 43 where said mesh generator is TRUEGRID. 

48. A system as defined in claim 43 where said stress/strain/deformation analyzer is 
DYNA3D. 
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49. A system as defined in claim 43 where said stress/strain/deformation analyzer is 
NIKE3D. 

50. A system as defined in claim 48 where said DYNA3D is modified to accommodate a 
strain energy density of the form: 

5 W -a (7-3) + a (7-3) + a (7-3) 2 -f-a (7-3)(7-3) + a (I -3) 2 + 

10 V 1 7 01 V 2 7 20 V 1 7 U V 1 7V 2 7 02 V 2 7 

a (7 -3) + a (7 -3) \l -3) + a (7 -3)(7 -3) 2 + a (7-3)' + 

30 V 1 7 21 V 1 7 V 2 7 12 V 1 7V 2 7 03 v 2 7 

M(/-l) 2 

with K = 2(& + a )/(l-2v) 

v 10 01 7 v 7 

1 0 where 

a., are material parameters; 

v is Poisson's ratio; 

is the bulk modulus given as a function of Poisson's ratio; and 
J , 1 9 and 7 3 are the first, second, and third invariants of the right Cauchy-Green 

1 5 strain tensor, respectively. 

51. A system as defined in claim 49 where said NIKE3D is modified to accommodate a strain 
energy density of the form: 

W =a (7 -3) + a (7 -3) + a (7-3) 2 + a (7-3)(7-3) + a (I-3) 2 + 

10 v 1 7 01 v 2 7 20 v 1 7 ll v 1 7V 2 7 02 v 2 7 

a (7 -3) + a (7 -3) 2 (7 -3) + a (7 -3)(7 -3) 2 + a (7-3)' + 

30 v 1 7 21 v 1 7 v 2 7 12 v 1 7V 2 7 03 v 2 7 

20 l AK{I-\f 

with A>2(a +a )/(l-2v) 
v io or v 7 

where 
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a . are material parameters; 



v is Poisson's ratio; 



K is the bulk modulus given as a function of Poisson's ratio; and 



I , I , and I are the first, second, and third invariants of the right Cauchy-Green 



53. 

10 54. 
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tensor, respectively. 

A system as defined in claim 43 further comprising visualization tool that receives said 
stresses and strains on said medical device and anatomical feature and displays said 
stresses and strains on said medical device via visual representation. 
A system as defined in claim 52 where said visualization tool is GRIZ. 
A computer method for analyzing a medical device comprising: 

a) acquiring three-dimensional volumetric data of at least one anatomical feature; 

b) generating a geometric model of said three-dimensional volumetric data; 

c) receiving data representing a geometric model of a candidate medical device design; 

d) receiving said geometric model of said three-dimensional volumetric data; 

e) generating a mesh incorporating both said geometric model of said anatomical feature 
and said geometric model of said candidate medical device design; 

f) receiving material properties of said mesh; 

g) receiving load data of said mesh; and 

h) simulating stresses, strains, and deformation imposed on said candidate medical device 
design by said load data. 

A method as defined in claim 54 further comprising the step of simulating stresses, 
strains, and deformations to a point of failure of said candidate medical device design. 
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56. A method as defined in claim 54 where said three-dimensional volumetric data are 
acquired via CT scan. 

57. A method as defined in claim 54 where said three-dimensional volumetric data are 
acquired via MRL 

5 58. A method as defined in claim 54 where said geometric model of a medical device is for 
an endovascular prosthesis. 

59. A method as defined in claim 58 where said endovascular prosthesis is a transluminally 
placed endovascular graft. 

60. A method as defined in claim 59 where said endovascaular prosthesis is a cardiovascular 
1 0 stent device. 

61 . A method as defined in claim 54 where said geometric model for three-dimensional 
volumetric data is generated by a MIMICS software application. 

62. A method as defined in claim 54 where said mesh is generated by TRUEGRID. 

63. A method as defined in claim 54 where said stresses, strains, and deformations are 
1 5 simulated by a DYNA3D software application. 

64. A method as defined in claim 54 where said stresses, strains, and deformations are 
simulated by a NIKE3D software application. 

65. A method as defined in claim 63 where said DYNA3D is modified to accommodate a 
strain energy density of the form: 

20 W =a (7-3) + a (7-3) + a (7-3) 2 + a (7-3)(7-3) + a (I -3) 2 + 

10 V 1 J 01 V 2 J 20 V 1 J 11 V 1 /V 2 J 02 V 2 y 

a (7 -3) + a (7 -3) \l -3) + a (7 -3)(7 -3) 2 + a (7-3) 3 + 

30 V 1 * 21 V 1 7 V 2 y 12 V 1 yV 2 7 03 v 2 y 
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with AT=2(a +a )/(l-2v) 

v 10 or v 7 

where 

a are material parameters; 

5 v is Poisson's ratio; 

K is the bulk modulus given as a function of Poisson's ratio; and 

I , I , and / are the first, second, and third invariants of the right Cauchy-Green 

strain tensor, respectively. 

66. A method as defined in claim 64 where said NIKE3D is modified to accommodate a 
10 strain energy density of the form: 

W =a (7-3) + a (7-3) + a (7-3) 2 + a (7-3)(7-3) + a (I -3) 2 + 

10 V 1 7 01 v 2 7 20 y 1 7 11 v 1 7V 2 7 02 V 2 7 

a (7 -3) + a (7-3) 2 (7-3) + a (7 -3)(7 -3) 2 + a (7-3) 3 + 

30 V 1 7 21 V 1 7 V 2 7 12 V 1 7V 2 7 03 V 2 7 

W-l) 2 



with 7T = 2(a +a )/(l-2v) 

V 10 01 7 v 7 



where 



a^ are material parameters; 

v is Poisson's ratio; 

is the bulk modulus given as a function of Poisson's ratio; and 
20 1 9 1 2 , and 1^ are the first, second, and third invariants of the right Cauchy-Green 

strain tensor, respectively. 

67. A method as defined in claim 54 where said stress/strain/deformation analysis is done 
using a non- linear finite element analysis tool. 
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A method as defined in claim 54 further comprising the step of receiving results of said 
stress, strain, and deformation analysis into a visualization tool and where said 
visualization tool visually presents the strains, stresses, and deformations on said medical 
device. 

A method as defined in claim 68 where said visualization means is GRIZ. 
A method for analyzing a medical device comprising: 

a) acquiring three-dimensional volumetric data of at least one anatomical feature of a 
particular individual; 

b) generating a geometric model of said three-dimensional volumetric data; 

c) receiving a geometric model of a candidate medical device; 

d) receiving said geometric model of said three-dimensional volumetric data; 

e) generating a mesh incorporating both said geometric model of said anatomical feature 
and geometric model of said candidate medical device; 

f) receiving material properties of said mesh; 

g) receiving load of said mesh; and 

h) simulating dynamic or quasi-static stresses, strains, and deformations imposed on 
medical device. 

A method as defined in claim 70 further comprising the step of simulating stresses, 
strains, and deformations to point of failure of said medical device. 
A method as defined in claim 70 where said three-dimensional volumetric data are 
acquired via CT scan. 

A method as defined in claim 70 where said three-dimensional volumetric data are 
acquired via MRI. 
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74. A method as defined in claim 70 where said geometric model of a medical device is for 

an endovascular prosthesis. 
75 A method as defined in claim 74 where said endovascular prosthesis is a transluminally 

placed endovascular graft. 
5 76. A method as defined in claim 74 where said endovascular prosthesis is a cariovascaula 

stent device. 

77. A method as defined in claim 70 where said generating geometric means for three- 
dimensional volumetric data is MIMICS. 

78. A method as defined in claim 70 where said mesh generating means is TRUEGRID. 

10 79. A method as defined in claim 70 where said stress/strain/deformation simulating means is 



80. A method as defined in claim 70 where said stress/strain/deformation simulating means is 



NIKE3D. 

81. A method as defined in claim 79 where said DYNA3D is modified to accommodate a 



DYNA3D. 



15 



strain energy density of the form: 



W 



= a (7-3) + a (7-3) + a (7-3) 2 + a (7-3)(7-3) + a (I-3) 2 + 

10 v 1 7 01 v 2 y 20 v 1 J 11 v 1 /x 2 J 02 v 2 ' 

a (7-3) + a (/ -3) \l -3) + a (7 -3)(7 -3) 2 + a (7-3) 3 + 

30 v 1 ' 21 V 1 ' V 2 ' 12 v 1 /v 2 J 03 v 2 ' 



with 




20 



where 



a., are material parameters; 



v is Poisson's ratio; 
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10 



K is the bulk modulus given as a function of Poisson's ratio; and 

7 , I , and / are the first, second, and third invariants of the right Cauchy-Green 

strain tensor, respectively. 

82. A method as defined in claim 80 where said NIKE3D is modified to accommodate a 
strain energy density of the form: 

W =a (7-3) + a (7-3) + a (7-3) 2 + a (7-3)(7-3) + a (I-3) 2 + 

10 v 1 7 01 v 2 7 20 v 1 7 ll v 1 7V 2 7 02 v 2 7 

a (7 -3) + a (7-3) 2 (7-3) + a (7 -3)(7 -3) 2 + a (7-3) 3 + 

30 v 1 7 21 v 1 7 v 2 7 12 v 1 7V 2 7 03 v 2 7 

with 7T=2(a +a )/(l -2v) 
v 10 or v 7 

where 



a are material parameters; 
vis Poisson's ratio; 

K is the bulk modulus given as a function of Poisson's ratio; and 
1 5 1 , 7 2? and 7^ are the first, second, and third invariants of the right Cauchy-Green 

strain tensor, respectively. 

83. A method as defined in claim 70 where said stress/strain/deformation analysis is done 
using a non-linear finite element analysis tool. 

84. A method as defined in claim 70 further comprising the step of receiving results of said 
20 stress and strain analysis into a visualization tool and where said visualization tool 

visually presents the strains and stresses on said medical device. 

85. A method as defined in claim 84 where said visualization means is GRIZ. 
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86. A computer method for analyzing a medical device comprising: 

a) receiving data representing an in vitro model and a geometric model of a candidate 
medical device design; 

e) generating a mesh incorporating both said geometric model of said in vitro model and 
5 geometric model of said candidate medical device design; 

f) receiving material properties of said mesh; 

g) receiving load data of said mesh; and 

h) simulating stresses, strains, and deformations imposed on said medical device by said 
load data. 

10 87. A method as defined in claim 86 farther comprising the step of simulating stresses and 
strains to point of failure of said medical device. 

88. A method as defined in claim 86 where said geometric model of a medical device is for 
an endovascular prosthesis. 

89. A method as defined in claim 86 where said endovascular prosthesis is a transluminally 
1 5 placed endovascular graft. 

90. A method as defined in claim 88 where said endovascular prosthesis is a cardiovascular 
stent device. 

91. A method as defined in claim 86 where said mesh generating means is TRUEGRID. 

92. A method as defined in claim 86 where said stress/strain/deformation simulating means is 
20 DYNA3D. 

93. A method as defined in claim 86 where said stress/strain/deformation simulating means is 
NIKE3D. 
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94. A method as defined in claim 92 where said DYNA3D is modified to accommodate a 
strain energy density of the form: 

W -a (7-3) + a (7-3) + a (/-3) 2 + a (7-3)(7~3) + a (I-3) 2 + 

10 I 01 v 2 7 20 v 1 7 11 v 1 7V 2 7 02 v 2 7 

a (7 -3) + a (7-3) 2 (7-3) + a (7 -3)(7 -3) 2 + a (7-3) 3 + 

30 V 1 7 21 V 1 7 V 2 7 12 V 1 7V 2 7 03 V 2 7 

with K = 2U +a )/(l-2v), 

v 10 or v 75 



where 

a. are material parameters; 

10 v is Poisson's ratio; 

K is the bulk modulus given as a function of Poisson's ratio; and 

7 , 7 2? and 7 3 are the first, second, and third invariants of the right Cauchy-Green 

strain tensor, respectively. 

95. A method as defined in claim 93 where said NIKE3D is modified to accommodate a 
1 5 strain energy density of the form: 

W =a (7-3) + a (7 -3) + a (7-3) 2 + a (7-3)(7-3) + a (I -3) 2 + 

10 v 1 7 01 v 2 7 20 v 1 7 ll v 1 7V 2 7 02 V 2 7 

a o (/-3) + a 2i (/-3) 2 (/ -3) + a 2 (7 -3)(7 -3) 2 + a^/ -3) 3 + 
with K = 2(a o + a m ) / (1 - 2v), 

20 

where 
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a are material parameters; 

y 

v is Poisson's ratio; 

Kis the bulk modulus given as a function of Poisson's ratio; and 

/ I and I are the first second, and third invariants of the right Cauchy-Green 

l' 2 ? 3 

strain tensor, respectively. 

96. A method as defined in claim 86 where said stress/strain/deformation analysis is done 
using a non-linear finite element analysis tool. 

97. A method as defined in claim 86 further comprising the step of receiving results of said 
stress, strain, and deformation analysis into a visualization tool and where said 
visualization tool visually presents the strains and stresses on said medical device. 

98. A method as defined in claim 97 where said visualization means is GRIZ. 

99. A method for analyzing a medical device comprising: 

a) receiving a geometric model of a candidate medical device design; 

b) generating a mesh of said candidate medical device design; 

c) receiving material properties of said mesh; 

d) receiving load of said mesh; and 

e) simulating stresses, strains, and deformations imposed on said medical device. 

100. A method as defined in claim 99 further comprising the step of simulating stresses and 
strains to point of failure of said medical device. 

101. A method as defined in claim 99 where said geometric model of a medical device is for 
an endovascular prosthesis. 
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102. 

103. 

5 104. 
105. 

106. 

10 107. 



A method as defined in claim 99 where said endovascular prosthesis is a transluminally 
placed endovascular graft. 

A method as defined in claim 101 where said endovascular prosthesis is a cardiovascular 
stent device. 

A method as defined in claim 99 where said mesh generating means is TRUEGRID. 

A method as defined in claim 99 where said stress/strain/deformation simulating means is 



A method as defined in claim 99 where said stress/strain/deformation simulating means is 
NIKE3D. 

A method as defined in claim 105 where said DYNA3D is modified to accommodate a 
strain energy density of the form: 



where 

a are material parameters; 
v is Poisson's ratio; 

K is the bulk modulus given as a function of Poisson's ratio; and 

J , 7, and / are the first, second, and third invariants of the right Cauchy-Green 

tensor, respectively. 



W 



-a (J -3) + a (7-3) + a (7-3) 2 + a <7-3)<7-3) + a (I-3) 2 + 

10 v 1 7 01 v 2 7 20 v l 7 ll v 1 7V 2 7 02 v 2 7 

a (7-3) + a (I -3) \l -3) + a (I -3)(/ -3) 2 + a (7-3) 3 + 

30 v 1 7 21 v 1 7 v 2 7 12 v 1 7V 2 7 03 v 2 7 

y^(/-i) 2 



with 
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108. A method as defined in claim 106 where said NIKE3D is modified to accommodate a 
strain energy density of the form: 

W =a (7-3) + a (7 -3) + a (7-3) 2 + a (7-3)(7-3) + a (I-3) 2 + 

10 v 1 01 v 2 J 20 v 1 7 ll v 1 ' v 2 7 02 v 2 7 

a (7 -3) + a (7 -3) \l -3) + a (7 -3)(7 -3) 2 + a (7-3) 3 + 

30 V 1 7 21 V 1 7 V 2 7 12 V 1 7V 2 7 03 y 2 7 

v^(/-i) 2 

with 7T=2(a +a )/(l-2v), 

v to or v 7? 

where 

a are material parameters; 
v is Poisson's ratio; 

K is the bulk modulus given as a function of Poisson's ratio; and 

1^, 1^ and 7^ are the first, second, and third invariants of the right Cauchy-Green 

strain tensor, respectively. 

109. A method as defined in claim 99 where said stress/strain/deformation analysis is done 
using a non-linear finite element analysis tool. 

110. A method as defined in claim 99 further comprising the step of receiving results of said 
stress, strain, and deformation analysis into a visualization tool and where said 
visualization tool visually presents the strains and stresses on said medical device. 

111. A method as defined in claim 110 where said visualization means is GRIZ. 
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Abstract 

A system and method of developing better-designed medical devices, particularly 
5 cardiovascular stents and endovascular grafts. The system comprises a geometry generator, a 
mesh generator, a stress/strain/deformation analyzer, and a visualization tool. In one 
embodiment, the geometry generator receives three-dimensional volumetric data of an 
anatomical feature and generates a geometric model. The mesh generator then receives such 
geometric model of an anatomical feature or an in vitro model and a geometric model of a 

10 candidate medical device. In another embodiment, the mesh generator only receives a geometric 
model of the candidate medical device. Using the geometric model(s) received, the mesh 
generator creates or generates a mesh or a finite element model. The stress/strain/deformation 
analyzer then receives the mesh, and the material models and loads of that mesh. Using analysis, 
preferably non-linear analysis, the stress/strain/deformation analyzer determines the predicted 

1 5 stresses, strains, and deformations on the candidate medical device. Such stresses, strains, and 
deformations may optionally be simulated visually using a visualization tool. 
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FIG. 5A 

Line Command 

1 c *** Slotted Tube Integrated Stent Design Simulation: istent.run **** 

2 c 

3 c parameter settings 

4 c 

5 c .... inike=l => make nike file; inike-0 => make dyna file 

6 c .... imodel = 0 ==> full 3 segment model with interconnects 

7 c = 1 => 3 -crown segment only 

8 c = 2 => 6-crown segment only 

9 c = 3 => 12-crown segment only 

10 c .... isym = 0 => full 360 deg model 

11 c = 1 — > symmetric model 

12 c .... isim_mode: type of simulation 

13 c = 1 : => radial force to R_f= X% RJO, restoring stress mat'l 

14 c =2: => flat plate force, restoring stress mat! 

15 c = 3 : => predelivery compression, loading stress mat'l 

16 c = 4: => initial expansion 

17 c = 5: => frequency analysis 

18 c .... refine = X => add X elements via mseq in each direction 

19 c of the cross section 

20 c 

21 parameter inike 1 ; 

22 parameter imodel 0 ; 

23 parameter isym 0 ; 

24 parameter isim_mode 4 ; 

25 parameter refine 2 ; 

26 c 

27 para Tighten [0.9]; c helps 'tighten 1 or stiffen spline 

28 c range (0.5,1) (probably should not change) 

29 c 

30 c parameter settings 

31 c 

32 c - design parameters ===== 

33 c 

34 c Note: Adjust specified OD for each segment considering the wall 

35 c thickness for that segment so that ID's match in a consistent 

36 c way for the tube blank from which they were cut. 

37 c 

38 c Upper segment — 3 crowns 

39 c Middle segment — 6 crowns 




Line Command 

40 c Lower segment — 12 crowns (conical) 

41 c 

42 c Parameters for 3 -crown segment 

43 c 

44 para 

45 RCyl3 [.5*2/25.4] 

46 dCIA3 [-.00] c delta of center of inner arc for 3 crown segment (-:0) 

47 dCOA3 [0] c delta of center of outer arc for 3 crown segment (0:+) 

48 CW3 [.007] c Circumferential width of segments for 3 crowns 

49 RW3 [,005] c Radial width for 3 crowns 

50 NRA3 [.0095] c normal radius of smaller cylinders (arcs) 

51 c for 3 crowns 

52 Ht3 [0.224] c distance from center of upper arcs 

53 c to center of lower arcs for 3 crowns 

54 NLegEB [12]; c number of elements along the leg 
55 

56 c 

57 c Parameters for 6-crown segment 

58 c 

59 para 

60 RCyl6 [,5*2/25,4] c outside radius for 6 crown segment 

6 1 dCIA6 [0] c delta of center of inner (smaller) arc for 6 crown 

segment(-:0) 

62 dCOA6 [0.002] c delta of center of outer (larger) arc for 6 crown 

segment (0:+) 

63 CW6 [.009] c Circumferential width of segments for 6 crowns 

64 RW6 [.009] c Radial width for 6 crowns 

65 NRA6 [.0105] c normal radius of smaller cylinders (arcs) 

66 c for 6 crowns 

67 Ht6 [.115] c distance from center of upper arcs 

68 c to center of lower arcs for 6 crowns 

69 NLegE16 [12]; c number of elements along the leg 
70 

71 c 

72 c Parameters for 12-crown segment 

73 c 

74 para 

75 dCIA12 [0] c delta of center of inner arc for 12 crown segment (-:0) 
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FIG. 5C 

Line Command 



76 dCOA12 [0] c delta of center of outer arc for 12 crown segment 

(0:+) 

77 CW12 [.005] c Circumferential width of segments for 12 crowns 

78 RW12 [.008] c Radial width for 12 crowns 

79 NRA12 [.006] c normal radius of smaller cylinders (arcs) 

80 c for 12 crowns 

81 Htl2 [.050] c distance from center of upper arcs 

82 c to center of lower arcs for 12 crowns 

83 c (measured along the leg, not necessarily in 

84 c the z direction) 

85 c first outside radius for 12 crown segment (near other segments) 

86 RCyll2_l [.5*2/25.4 - (.016-%RW12)] 

87 c second outside radius for 12 crown segment (bottom) 

88 RCY112_2 [.5* 1 .4/25.4 - (.016-%RW12)] 

89 c 

90 NLegE112 [10]; c number of elements along the leg 



91 

92 c 

93 c Interconnects 

94 c 
95 

96 c 

97 c Upper interconnects 

98 c 

99 paraHRJp [.02] c height of interconnect 



100 FRUp [.005] c fillet radius for blend 

101 ICWUp [.006] c circumferential width 

102 IRWUp3 [.005] c radial width at 3-crown end 

103 IRWUp6 [.006]; c radial width at 6-crown end 



104 

105 c 

106 c S-interconnects 

107 c 

108 para SIVer [.01] c vertical distance between upper or lower arc centers 



109 c also the distance from the vertical mid-line to 

110 c the first arc center 

111 SIHor [.010] c horizontal distance between upper two or 

112 c lower two arc centers 

113 Sir [.004] care radius 



Line Command 

114 SIrO [%SIr+%ICWUp/2] c outer radius 

1 1 5 SIrl [%SIr-%ICWUp/2]; c inner radius 
116 

117 c 

118 c Lower interconnects 

119 c 

120 paraHILr [.031] c height of interconnect 

121 FRLr [.010] c fillet radius for blend 

122 ICWLr [.007] c circumferential width 

123 IRWLr6 [.005] c radial width at 6-crown end 

124 IRWLrl2 [.005]; c radial width at 12-crown end 
125 

126 c 
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127 c .... = design parameters = 

128 c 

129 c .... set cylinder ID & OD for compression 

130 c 

13 1 if (%isim_mode.le.3) then 

132 parameter ricompcyl 

[l.l*max(%RCyl3 5 %RCyI6 5 %RCyll2_l 5 %RCyll2_2)] ; 

133 parameter rocompcyl 

[1.4*max(%RCyl3 9 %RCyl6,%RCyll2 - l 9 %RCyI12_2)] ; 

134 c 

135 c .... set cylinder ID & OD for expansion 

136 c 

137 elseif (%isim_mode.eq.4) then 

138 parameter rocompcyl 

[0.95*(min(%RCyl3,%RCyI6 5 %RCyll2_l 5 %RCyll2_2)-%RW6)] ; 

139 parameter ricompcyl 

[0.7*(min(%RCyi3,%RCyl6 5 %RCyll2_l,%RCyll2 - 2)-%RW6)] ; 

140 endif 

141 c 

142 c 

143 c Materials assignments 

144 c 

145 parameter matstl2 3 ; 

146 parameter matst6 4 ; 

147 parameter matst3 5 ; 



y... 



Line Command 

148 parameter matil26 6 ; 

149 parameter mati63 7; 

150 c 

151 c 

152 if (%isim_mode.eq. 1 ) then 

153 echo *** Radial Force Simulation *** 

154 elseif (%isim_mode.eq.2) then 

155 echo * * * Flat Plate Force Simulation * * * 

156 elseif (%isim_mode.eq.3) then 

157 echo *** Predelivery Compression Simulation *** 

158 elseif (%isim_mode.eq.4) then 

159 echo *** Initial Expansion Simulation *** 

160 elseif (%isim_mode.eq.5) then 

161 echo * * * Natural Frequency Analysis * * * 

162 else 

163 echo ! ! ! ERROR: illegal isim_mode ! ! ! 

164 interrupt 

165 endif 

166 c 

167 c analysis options 

168 title stent initial expansion simulation 

169 c 

1 70 c * * * DYNA3D Analysis Options * * * 

171 c 

172 if (%inike.eq.O) then 

1 73 echo Making DYNA3D input file 

174 dyna3d 

175 dynaopts 

176 term 5.0e-5 

177 plti l.e-6 

178 prti 5.0e-6 

179 c 

180 c .... DR options 

181 c 

182 itrx500 

183 tolrx 1.0e-2 

184 drdb 

185 c 

186 c .... thermal effects option - temp from load curve 1 
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FIG 5F 

Line Command x XVJ * ^ A 

187 c 

188 teo 1 

189 c 

190 tssfO.O 

191 c 

192 c print initial time step size 

193 c 

194 c prtflg 1 

195 c 

196 c .... turn off (0) or on (1) SAND database flag 

197 c 

198 edsdfO 

199 c 

200 nrest 90000 

201 nrunr 95000 ; 

202 c 

203 c .... DYNA3D discrete nodes impacting surface - stent to cyl 

204 c * one side (180 deg) * 

205 c 

206 sid 1 dni 

207 c sfif 

208 c mfif 

209 pnlts 1.0e-0 

210 pnltm 1.0e-0 

211 ; 

212 c 

213 c .... DYNA3D discrete nodes impacting surface - stent to cyl 

214 c * opposite side * 

215 c 

216 c sid 2 dni 

217 c sfif 

218 cmfif 

219 c pnlts 1.0e-4 

220 c pnltm 1.0e-4 

221 c ; 

222 c 

223 c .... end DYNA3D commands 

224 c 

225 endif 
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226 c 

227 c 

228 c * * * NIKE3D Analysis Options * * * 

229 c 

230 if (%inike.eq. 1) then 

23 1 echo Making NIKE3D input file . . . 

232 nike3d 

233 nikeopts 

234 nstep 5 

235 delt0.2 

236 anal stat 

237 c 

238 c .... step tol of le-8 seems OK for predel compression 

239 c 

240 if (%isim_mode.eq. 1 ,or.%isim_mode.eq.2) then 

241 dctol-1.0e-8 

242 elseif (%isim_mode.eq.3) then 

243 dctol-1.0e-6 

244 endif 

245 c 

246 c .... max iterations per stiffness reform 

247 c 

248 nibsr 20 

249 c 

250 c .... max stiffness reforms per step 

251 c 

252 msrf 20 ; 

253 c 

254 c .... temperatures follow load curve 1 

255 c ** manually add tref=1.0 on matl 2 control card cols 26-35 * 

256 c 

257 teo 1 

258 if (%isim_mode.eq. 1 .or.%isim_mode.eq.2) then 

259 iprt 1 

260 elseif (%isim_mode.eq.3.or.%isim_mode.eq.4) then 

261 iprt 25 

262 endif 

263 iplt 1 

264 nsbrr 1 



Line Command 

Itl *f core 1 FIG. 5H 

266 bfgscore 

267 bwmo new 

268 echo Bandwidth minimization ACTIVATED with "NEW" option 

269 c 

270 c element constitutive data incore 

271 c 

272 bfor 10 

273 sfor 10 

274 bef 11 

275 c 

276 c .... linear solver 

277 c 

278 lsolver fissle 

279 c 

280 c .... solid element stent contact surface 

281 c 

282 sid 1 sv 

283 c 

284 if(%isim_mode.eq.l)then 

285 c 

286 c 

287 pnlt 1.0e-5 

288 elseif (%isim_mode.eq.2) then 

289 pnlt 0.00001 

290 elseif (%isim_mode.eq.3) then 

291 c 

292 c .... essential to adjust penalty 

293 c 

294 pnlt 1.0e+4 

295 elseif (%isim_mode.eq.4) then 

296 pnlt 1.0e-5 

297 c iaug 1 ; 

298 endif 

299 ; 

300 c 

301 c .... slidesurface between interconnects and segments 

302 c 

303 sid 2 tied 






Line 
304 
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305 c 

306 c .... NIKE3D shell geometric stiffness (HL only) 

307 c 

308 segs 1 ; 

309 c 

310 c .... end NIKE3D section 

311 c 

312 endif 

313 c 

3 14 c .... symmetry planes 

315 c 

316 if (%isym.eq. 1) then 

317 c 

318 c... Symmetric Model 

319 c theta=-60 and +60 symmetry to remove rigid body modes 

320 c 

321 c plane 1 

322 c 0.0 0.0 0.0 

323 c [-sin(60)] [-cos(60)] 0.0 

324 c 0.0005 symm ; 

325 c plane 2 

326 c 0.0 0.0 0.0 

327 c [-sin(60)] [cos(60)] 0.0 

328 c 0.0005 symm ; 

329 c 

330 else 

331 c 

332 c .... symmetry planes to remove rigid body modes for full model 

333 c 

334 plane 1 

335 0.0 0.0 0.0 

336 1.0 0.0 0.0 

337 .0005 symm ; 

338 plane 2 

339 0.0 0.0 0.0 

340 0.0 1.0 0.0 

341 .0005 symm ; 

342 c plane 3 



'V 



if 



Line Command 

343 c 0.0 0.0 0.0 FTP ST 

344 c 0.0 0.0 TBD r\\J. J J 

345 c .0005 symm ; 

346 endif 

347 c 

348 c 

349 if(%inike.eq.0)then 

350 c 

35 1 c .... Load Curves for DYNA3D **ADD DR FLAG TO INPUT FILE ** 

352 c 

353 if (%isim_mode.eq.l)tlien 

354 c 

355 c .... radial force 

356 c 

357 led 1 

358 O.000E+00 LOOOE-KK) 

359 7.500E-03 2.250E+04 

360 1 .000E-00 2.250E+04 ; 

361 c 1.000E-02 3.000E+04 

362 c 1.000E-00 3.000E+04 ; 

363 elseif (%isim_mode.eq.2) then 

364 c 

365 c .... flat plate compression, led 1 not used (dummy definition) 

366 c 

367 quit 

368 c 

369 elseif (%isim_mode.eq.3) then 

370 c 

371 c . . . . predelivery compression strain 

372 c 

373 led 1 

374 O.OOOE+OO 1.000E+00 

375 1.000E-02 2.008E+05 

376 1.000E-00 2.008E+05 ; 

377 endif 

378 c 

379 c .... load curve #2 only used for flat plate compression 

380 c 

381 led 2 
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Line Command 

382 0.000E+00 0.000E+00 FIG. 5K 

383 1.000E+00 0.000e-00 ; 

384 endif 

385 c 

3 86 if (%inike.eq. 1 ) then 

387 c 

388 c .... ******** Load Curves for NIKE3D ********** 

389 c 

390 if (%isim_mode.eq. 1) then 

391 c 

392 c .... radial force 

393 c 

394 led 1 

395 0.000E+00 1.000E+00 

396 1.000E+00 2.000E+03 ; 

397 elseif (%isim_mode.eq.2) then 

398 c 

399 c .... flat plate compression 

400 c 

401 led 1 

402 0.000E+00 1.000E+00 

403 1.000E+00 0.000E+00 ; 

404 elseif (%isim_mode.eq.3) then 

405 c 

406 c .... predelivery compression strain 

407 c 

408 led 1 

409 0.000E+00 1.000E+00 

410 1.000E+00 2.008E+03 ; 

411 elseif (%isim_mode.eq.4) then 

412 c 

413 c .... initial expansion strain 

414 c 

415 led 1 

416 c .... thermal load (activate TEO above) 

417 c 0.000E+00 1.000E+00 

418 c 1.000E+00-2.008E+04 ; 

419 c . . . . prescribed displacement 

420 0.000E+00 O.OOOE+OO 
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Line Command 

421 1.000E+00 L000E-02; uyp <-t 

422 endif riVJ - JLj 

423 c 

424 c stent parts 

425 c 

426 include irss.tg 

427 c 

428 c stent materials 

429 c 

430 if (%inike.eq. 1) then 

43 1 if (%isirnjmode.eq. 1 .or.%isim_mode.eq.2) then 

432 include istent.mats_nike_solid 

433 echo model for radial force/flat plate analysis 

434 elseif (%isim_mode.eq.3) then 

43 5 include istent jnats_compress_nike_solid 

436 echo model for predelivery compression strain 

437 elseif (%isiin_mode.eq.4) then 

438 include istentmate_compress_nike_solid 

439 echo model for initial expansion strain 

440 endif 

441 c 

442 elseif (%inike.eq.O) then 

443 if (%isim_mode.eq. l,or.%isim_mode.eq.2) then 

444 include istentmats_dyna_solid 

445 echo model for radial force/flat plate analysis 

446 elseif (%isim_mode.eq.3) then 

447 include istent.mats_ compress^ dyna_solid 

448 echo model for predelivery compression strain 

449 elseif (%isim_mode.eq.4) then 

450 include istentmats_compress_dyna_solid 

45 1 echo model for initial expansion strain 

452 endif 

453 endif 

454 c 

455 c .... cylindrical compression for radial force or predelivery compression 

456 c 

457 if (%isim_mode.eq. 1 .or.%isim_mode.eq.3 .or.%isim_mode.eq.4) then 

458 c 

459 if (%isym.eq.l) then 



Line Command 

460 include cylinder.parts sym 

461 else 

462 include cylinder .parts 

463 endif 

464 c 

465 if (%inike.eq. 1 ) then 

466 include cylinder.materialsnike 

467 elseif (%inike.eq.O) then 

468 include cylinder.materials_dyna 

469 endif 

470 endif 

471 c 

472 stp .01 

473 merge 
.474 c 
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1 c ********** TPEG Inflatable Proximal Seal Simulation ****************** 

2 c (seaLrun) 

3 c March, 1999 JrlLr. OA 

4 c 

5 c parameter settings 

6 c 

7 c .... analytical model aorta geometric parameters 

8 c (distortion is 4-lobe) 

9 c 

1 0 parameter r_aorta [10.0/25.4] ; 

1 1 parameter thk_aorta [1.0/25.4] ; 

12 parameter amp_plaque [0.0/25.4] ; 

13 c 

14 parameter ro_aorta [%r_aorta+%fhk_aorta] ; 

15 c 

16 c .... - TPEG Design Parameters - 

17 c 

18 parameter t_tpeg [10/25.4] ; 

19 parameter r_ps [3/25.4] ; 

20 parameter l_tpeg 2.0 ; 

21 parameter l_flap 0.25 ; 

22 c 

23 parameter graft_wall_thick [6*0.0013] 

24 parameter cuff_wall_thick [3*0.0013] 

25 parameter flap_wall_thick [6*0.0013] 

26 c 

27 c 

28 c .... Pressures and load curve assignments 

29 c 

30 parameter P_hemo 2.32 ; 

3 1 parameter P_cuff 3.0 ; 

32 c 

33 parameter Ic_hemo 1 ; 

34 parameter lc_proxcuff 3 ; 

35 c 

36 c .... TPEG folding simulation parameters 

37 c 

38 parameter vel_fold 20.0 ; 

39 parameter t_fold [0.25/%vel_fold] ; 

40 parameter t init 0.0e-3 ; 

41 c 

42 c 
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43 c analysis options 

44 title sc6.i Seal CT-Solid r_t=lOmm rjps=3mm P_cuff=3.0 990428 FIG. 6B 

45 c 

46 c *** DYNA3D Analysis Options *** 

47 c 

48 dyna3d 

49 dynaopts 

50 term 6.5e-2 

51 plti5.e-4 

52 prti 2.5e-2 

53 c 

54 c .... DR options 

55 c 

56 itrx 500 
□ 57 c 

i -■' 58 c .... increase DR tol to prevent convergence after compression before expansion 

r\ 59 c 

J 60 ctolrxl.0e-6 

v| 61 tolrxl.0e-12 

^ 62 drdb 

/ ! 63 c 

ri 64 tssf0.9 

3 65 c 

J3 66 c .... turn off (0) or on (1) SAND database flag 

£J 67 c 

68 edsdfO 

" 69 c 

70 nrest 90000 

71 nrunr 5000 ; 

72 c 

73 c .... symmetry planes on xz and yz planes 

74 c 

75 plane 1 

76 0.0 0.0 0.0 

77 1.0 0.0 0.0 0.001 symm; 

78 plane 2 

79 0.0 0.0 0.0 

80 0.01.0 0.0 0.001 symm; 

81 c 

82 c .... DYNA3D slidesurface: +x folder cylinder 

83 c 

84 sid 1 sv 




85 pnlts 1.0 

86 pnltm 1.0 FIG. 6C 

87 pen 

88 ; 

89 c 

90 c .... DYNA3D slidesurface: -x folder cylinder 

91 c 

92 sid 2 sv 

93 pnlts 1.0 

94 pnltm 1.0 

95 pen 

96 ; 

97 c 

98 c 

99 c .... DYNA3D slidesurface: +y folder cylinder 

100 c 

101 sid 3 sv 

102 pnlts 1.0 

103 pnltm 1.0 

104 pen 

105 ; 

106 c 

107 c .... DYNA3D slidesurface: -y folder cylinder 

108 c 

109 sid 4 sv 

110 pnlts 1.0 

111 pnltm 1.0 

112 pen 

113 ; 

114 c 

115 c .... DYNA3D tpeg to aorta (aorta is master) 

116 c 

117 sid5sv 

118 c 

119 c .... solid element aorta 

120 c 

121 pnlts 0.1 

122 pnltm 0.1 

123 c 

124 c .... shell element aorta 

125 c 

126 c pnlts 1.0 




127 cpnltml.O t^t/^ r ta 

128 pen FIG. 6D 

129 ; 

130 c 

131 c .... load curve: hemodynamics **** ADD DR FLAG TO INPUT FILE **** 

132 c 

133 led 1 

134 0.000E+00 0.000E+00 

135 [%t_init+2*%t_fold+1.0e-3] 0.000e+00 

136 [%t_init+2*%t_fold+2.0e-3] %P_hemo 

137 1.000E+00 %P_hemo ; 

138 c 

139 c..„ load curve: channel !! NOT USED !!**** ADD DR FLAG TO INPUT FILE **** 

140 c 

141 led 2 

142 0.000E+O0 O.0O0E+00 

143 [%t_init+2*%t_fold+1.0e-3] 0.000e+00 

144 [%t_init+2*%t_fold+2.0e-3] 0.000e-00 

145 1.000E+00 0.000e-00; 

146 c 

147 c .... load curve: proximal cuff **** ADD DR FLAG TO INPUT FILE **** 

148 c 

149 led 3 

150 0.O00E+00 0.000E+00 

151 [%t_init+2*%t_fold+1.0e-3] 0.000e+00 

152 [%t_init+2*%t_fold+2.0e-3] %P_cuff 

153 1.000E+00 %P_cuff; 

154 c 

1 55 c .... load curve for +x folder cylinder motion/velocity 

156 c 

157 led 4 

158 0.000E+00 0.000E+00 

159 %t_init 0.000E+00 

160 [%t_init+1.0E-04] [-%vel_fold] 

161 [%t_init+%t_fold] [-%vel_fold] 

162 [%t_init+%t_fold+1.0e-3] 0.000E+00 

163 [%t_init+2*%t_fold+1.0e-3] 0.000e+00 

164 [%t_init+2*%t_fold+2.0e-3] [2.0*%vel_fold] 

165 [%t_init+3*%t_fold+2.0e-3] [2.0*%vel_fold] 

166 [%t_init+3*%t_fold+3.0e-3] 0.000e+00 

167 1.000E+00 0.000E+00 ; 

168 c 
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169 c .... load curve for -x folder cylinder motion 



176 [%t_init+%t_fold+1.0e-3] 0.000E+00 

177 [%t_init+2*%t_fold+1.0e-3] 0.000e+00 

178 [%t_init+2*%t_fold+2.0e-3] [-2.0*%vel_fold] 

179 [%t_init+3*%t_fold+2.0e-3] [-2.0*%vel_fold] 

180 [%t_init+3*%t_fold+3.0e-3] 0.000e+00 

181 1.000E+00 0.000E+00 ; 

182 c 

1 83 c .... load curve for +y folder cylinder motion 

184 c 

185 led 6 

186 0.000E+00 O.OOOE+00 

187 %t_init 0.000E+00 

188 [%t_init+1.000E-04] [-%vel_fold] 

189 [%t_init+%t_fold] [-%vel_fold] 

190 [%t_init+%t_fold+1.0e-3] O.OOOE+00 

191 [%t_init+2*%t_fold+1.0e-3] 0.000e+00 

192 [%t_init+2*%t_fold+2.0e-3] [2.0*%vel_fold] 

193 [%t_init+3*%t_fold+2.0e-3] [2.0*%vel_fold] 

194 [%t_init+3*%t_fold+3.0e-3] 0.000e+00 

195 1.000E+00 O.OOOE+00 ; 

196 c 

197 c .... load curve for -y folder cylinder velocity 

198 c 

199 led 7 

200 0.000E+00 0.000E+00 

201 %t_init 0.000E+00 

202 [%t_init+1.000E-04] [ %vel_fold] 

203 [%t_init+%t_fold] [ %vel_fold] 

204 [%t_init+%t_fold+1.0e-3] 0.000E+00 

205 [%t_init+2*%t_fold+1.0e-3] 0.000e+00 

206 [%t_init+2*%t_fold+2.0e-3] [-2.0*%vel_fold] 

207 [%t_init+3*%t_fold+2.0e-3] [-2.0*%vel_fold] 

208 [%t_init+3*%t_fold+3.0e-3] 0.000e+00 

209 1.000E+00 0.O00E+00; 

210 c 



170 
171 
172 
173 
174 
175 



c 
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led 5 
0.000E+00 
%t_init 

[%t_init+1.000E-04] 
[%t_init+%t_fold] 



0.000E+00 
0.O00E+0O 

[ %vel_fold] 
[ %vel_fold] 
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211 
212 
213 



c 



parts and materials 



FIG. 6F 



c 



c 



214 c .... get CT-data meshed aorta; convert cm to inches 

215 c 

216 csca[l./2.54] 

217 include tpeg.part_ct_aorta3 

218 c 

219 csca 1.0 

220 c 

221 c .... option for analytical aorta model 

222 c 

223 c include tpeg.part_eq_aorta 

224 c 

225 include tpeg.part_cuffl 

226 include tpeg.part_folder2 

227 c 

228 include tpeg.materials_dyna 

229 c 

230 c .... use negative tols to prevent aorta nodes merging w/ folder cylinder 

231 c nodes if they coincidently become adjacent 

232 c 

233 c .... merge nodes within CT aorta part using rather loose tolerance 

234 c 

235 bptol 1 1 0.01 

236 bptol 1 3 -1.0 

237 bptol 14 -1.0 

238 bptol 1 5 -1.0 

239 bptol 16 -1.0 

240 tp .001 

241 c 



1 c FIG. 7A 

2 c tpeg.part_ct_aorta3 

3 c April 15, 1999 

4 c 

5 c Aortic Model for Inflatable TPEG Model 

6 c Derived from Patient CT Data 

7 c Outer surface constructed with 0.52 mm offset from inner 

8 c 

9 c .... this is an aortic mesh file which surrounds the neck of the 

10 c 3-D AAA reconstruction with solid elements. 

11 c 

12 c This file uses TrueGrid planes, oriented by eye using trial 

13 c and error graphically, to determine an orthonormal section. 

14 c Trick there is to adjust surface until walls of proximal neck section 

15 c are parallel to global z axis. Use rz to rotate screen to find values, 

16 c then use in surface transformation to position CT data for meshing. 

17 c 

18 c .... import IGES file containing surface data from CT scan 

19 c 

20 iges solidl .igs 1 1 mx -18.54 my -16.8 ry 24 rx 22 mz 4.8; 

21 c 

22 c .... inner surface 

23 c 

24 sd 17 sds 9 12; 

25 c 

26 c .... outer surface 

27 c 

28 sd 18 sds 15 16 ; 

29 c 

30 sd 201 plan 

31 0.0.1.5 

32 0 0 1 

33 sd 202 plan 

34 0. 0. 2.5 

35 0 0 1 

36 sd 203 plan 

37 0. 0. -2.3 

38 0 0 1 

39 sd 204 plan 

40 0. 0. 3.3 

41 0 0 1 

42 sd301 cyOOOOO 1 1.35 





43 
44 
45 
46 
47 
48 
49 
50 



sd 401 plan 
0. 0. 0. 
0. 1.0. 



c 



c .... adjust mz to position part at cuff on Z-axis; 
c cuff may be z=[2,2. 1 5] 
cylinder 
12; 
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51 12 3; 

52 12 3 4; 

53 c 

54 1.0 1.25 

55 0 180.0 360.0 

56 -2.3 1.5 2.5 3.3 
O 57 c 

58 mseqi2 

f: 59 mseqj29 29 

J 60 mseqk20 5 5 

sj 61 c 

V,, 62 c project top and bottom ends of aorta segment onto orthonormal planes 

63 c 

j\ 64 sfi;;-2; sd201 

5 65 sfi;;-3; sd202 

5 66 c 

Jp 67 c .... project top of upper neck segment onto orthonormal plane 

5 68 c 

u 69 sfi;;-4; sd204 

70 c 

71 c .... project bottom of lower neck segment onto orthonormal plane 

72 c after radially expanding bottom ring by delta-r=2.0 

73 mbi-1; ;-l;x2.0 

74 mbi-2; ;-l;x2.0 

75 sfi; ;-l; sd203 

76 c 

77 c .... project inner cylinder surface onto aorta luminal surface 



78 
79 
80 
81 
82 
83 
84 



c 



c .... project outer cylinder onto aorta outer wall surface 



c 



sfi-1; 1 3;2 3;sdl7 
sfi -1; 13; 3 4; sd 17 
sfi -1; 1 3; 1 2; sd 17 



c 




85 sfi -2; 1 3; 2 3; sd 18 

86 sfi -2; 1 3;3 4;sdl8 FIG. 7C 

87 sfi -2; 1 3; 1 2; sd 18 

88 c 

89 c .... project theta=0/360 seam onto a plane to facilitate merging 

90 c 

91 sfi 1 2;-l; ; sd401 

92 sfi 1 2; -3; ; sd 401 

93 c 

94 c 

95 c ... — slidesurface definition with TPEG body — 

96 c 

97 orpt + 0. 0. 3.0 

98 sii -1; 1 3; 3 4; 5 m 

99 c 

100 c .... +y hemicylinder is material 1 1 ; -y is mat 12 

101 c 

- 102 rati; 12; 24; 11 

i 103 mti;2 3;24; 12 

J 104 c 

105 c .... rigid material for aneurysm sac 

106 c 

! 107 mti ; 1 3; 1 2; 13 

I 108 c 

= 109 c .... Boundary Conditions 

* 1 10 c * fix proximal end only in z 

111 c 

112 bi;;-4; dzl ; 

113 c 

114 c .... adjust mz to position aorta at cuff on Z-axis; 

115 c cuff may be z=[2,2. 15] 

116 lctl 

117 mz [1.01*2.54] mx 0.7; ; 

118 lrep 1 ; 

119 endpart 

120 c 




1 c ****** Slotted Tube Integrated Stent Design Simulation ******* "CTp 1 O A 

2 c (istentrun) r HJT. O/V 

3 c Stent design analysis & CT- Anatomy simulation 

4 c 

5 c parameter settings 

6 c 

7 c .... inike=l => make nike file; inike=0 => make dyna file 

8 c .... imodel = 0 => full 3 segment model with interconnects 

9 c = 1 => 3 -crown segment only 

10 c = 2 => 6-crown segment only 

11 c = 3 => 12-crown segment only 

12 c .... isym = 0 => full 360 deg model 

13 c = 1 => symmetric model 

14 c isim_mode: type of simulation 

15 c = 1 : => radial force to R_f = 80% RjO, restoring stress mat 1 1 

16 c = 2: => flat plate force, restoring stress mat! 

17 c = 3 : => predelivery compression to 12 F, loading stress matl 

18 c = 4: => initial expansion 

19 c =5: => frequency analysis 

20 c =6: => anatomy deployment 

21 c .... refine = X => add X elements via mseq in each direction 

22 c of the cross section 

23 c 

24 c !!! warning - only 1st 8 characters of variable unique !!!! 

25 c 

26 parameter inike 1 ; 

27 parameter imodel 2 ; 

28 parameter isym 0 ; 

29 parameter isim_mode 6 ; 

30 parameter refine 1 ; 

31 c 

32 para Tighten [0.9] ; c helps 'tighten 1 or stiffen spline 

33 c range (0.5,1) (probably should not change) 

34 c 

35 c parameter settings 

36 c 

37 c .... • design parameters ============ 

38 c 

39 c Note: Adjust specified OD for each segment considering the wall thickness 

40 c for that segment so that ID's match in a consistent way for the 

41 c tube blank from which they were cut. 

42 c 

43 c Upper segment — 3 crowns 

44 c Middle segment ~ 6 crowns 

45 c Lower segment — 12 crowns (could be conical) 

46 c 

47 c Parameters for 3 -crown segment 

48 c 

49 para 



18/ 



50 RCyl3 [29*0.5/25.4] 

5 1 dCIA3 [-.00] c delta of center of inner arc for 3 crown segment (-:0) r lvJ. 

52 dCOA3 [0] c delta of center of outer arc for 3 crown segment (0:+) 

53 CW3 [.020] c Circumferential width of segments for 3 crowns 

54 RW3 [.018] c Radial width for 3 crowns 

55 NRA3 [.0195] c normal radius of smaller cylinders (arcs) 

56 c for 3 crowns 

57 Ht3 [1 .048] c distance from center of upper arcs 

58 c to center of lower arcs for 3 crowns 

59 NLegE13 [12]; c number of elements along the leg 

60 c 

61 c Parameters for 6-crown segment 

62 c 

63 para 

64 RCyl6 [29*0.5/25.4] c outside radius for 6 crown segment 

65 dCIA6 [0] c delta of center of inner (smaller) arc for 6 crown segment (-:0) 

66 dCOA6 [0.005] cdeltaof center of outer Garger) arc for 6 crown segment (0:+) 

67 CW6 [.020] c Circumferential width of segments for 6 crowns 

68 RW6 [.018] c Radial width for 6 crowns 

69 NRA6[.0195] c normal radius of smaller cylinders (arcs) 

70 c for 6 crowns 

71 Ht6 [.3 10] c distance from center of upper arcs 

72 c to center of lower arcs for 6 crowns 

73 NLegE16[12]; c number of elements along the leg 

74 c 

75 c Parameters for 12-crown segment 

76 c 

77 para 

78 dCIA12 [0] c delta of center of inner arc for 12 crown segment (-:0) 

79 dCOA12 [0] c delta of center of outer arc for 12 crown segment (0:+) 

80 CW12 [.008] c Circumferential width of segments for 12 crowns 

81 RW12 [.008] c Radial width for 12 crowns 

82 NRA1 2 [.006] c normal radius of smaller cylinders (arcs) 

83 c for 12 crowns 

84 Htl2 [-164] c distance from center of upper arcs 

85 c to center of lower arcs for 12 crowns 

86 c (measured along the leg, not necessarily in 

87 c the z direction) 

88 c first outside radius for 12 crown segment (near other segments) 

89 RCyll2_l [22*0.5/25.4] 

90 c second outside radius for 12 crown segment (bottom) 

91 RCY112_2 [20*0.5/25.4] 

92 c 

93 NLegEl 1 2 [ 1 0] ; c number of elements along the leg 

94 c 

95 c Interconnects 

96 c 

97 c Upper interconnects 

98 c 




99 para 

100 c HIUp [.10] c height of interconnect 

101 HIUp [.20] c height of interconnect 

102 FRUp [.016] c fillet radius for blend 

103 ICWUp [.010] c circumferential width 

104 IRWUp3 [.016] c radial width at 3-crown end 

105 IRWUp6 [,016]; c radial width at 6-crown end 

106 c 

107 c S-interconnects 

108 c 

109 para 

110 c SIVer [.03] c vertical distance between upper or lower arc centers 

1 1 1 SIVer [.06] c vertical distance between upper or lower arc centers 

112 c also the distance from the vertical mid-line to 

113 c the first arc center 

1 14 SIHor [.0125] c horizontal distance between upper two or 

115 c lower two arc centers 

116 Sir [.008] care radius 

117 SIrO [%SLS-°/oICWUp/2] c outer radius 

118 SM [%SIr-%ICWUp/2]; c inner radius 

119 c 

120 c Lower interconnects 

121 para 

122 c HELr [.071] c height of interconnect 

123 HILr [.142] c height of interconnect 

124 FRLr [.016] c fillet radius for blend 

125 ICWLr [.016] c circumferential width 

126 IRWLr6 [.005] c radial width at 6-crown end 

127 IRWLrl2 [.005]; c radial width at 12-crown end 

128 c 

129 c .... = design parameters 

130 c 

131 c .... set cylinder ID & OD for compression 

132 c 

133 if (%isim_mode.le.3 .or.%isim_mode.eq.6) then 

134 parameter ricompcyl [l.l*max(%RCyl3,%RCyl6 ? %RCyll2_l 5 %RCyll2_2)] ; 

135 parameter rocompcyl [1.4*max(%RCyl3 5 %RCyl6 ? %RCyll2_l ? %RCyll2_2)] ; 

136 c 

137 c .... set cylinder ID & OD for expansion 

138 c 

139 elseif (%isim_mode. eq.4) then 

140 parameter rocompcyl [0.95*(min( o /oRCyl3,%RCyl6 5 o /oRCyll2_l 5 %RCyll2„2)-%RW6)] ; 

141 parameter ricompcyl [0.7* (min(%RCyl3 5 %RCyl6 5 %RCyll2_l 5 %RCyll2_2)-%RW6)] ; 

142 endif 

143 c 

144 c Materials assignments 

145 c 

146 parameter matstl2 3 ; 

147 parameter matst6 4; 



FIG. 8C 




148 parameter matst3 5; 

1 49 parameter mati 1 26 6 ; F IG . 8 D 

150 parameter mati63 7; 

151 c 

152 if (%isim_mode.eq. 1) then 

153 echo *** Radial Force Simulation *** 

154 elseif (%isim_mode.eq.2) then 

155 echo *** Flat Plate Force Simulation *** 

156 elseif (%isim_mode.eq.3) then 

157 echo *** Predelivery Compression Simulation *** 

158 elseif (%isim_mode.eq.4) then 

159 echo *** Initial Expansion Simulation *** 

1 60 elseif (%isim_mode.eq.5) then 

161 echo * * * Natural Frequency Analysis * * * 

1 62 elseif (%isim_mode.eq.6) then 

163 echo *** Anatomy Deployment Simulation*** 

164 else 

165 echo !!! ERROR: illegal isimjnode !!! 

166 interrupt 

167 endif 

168 c 

1 69 c analysis options 

170 title human-size stent anatomy deployment 

171 c 

172 c *** DYNA3D Analysis Options *** 

173 c 

174 if (%inike.eq.O) then 

175 echo Making DYNA3D input file 

176 dyna3d 

177 dynaopts 

178 term2.0e-4 

179 pltil.e-4 

180 prti5.0e-6 

181 c 

182 c DR options 

183 c 

184 citrxSOO 

185 ctolre L0e-6 

186 cdrdb 

187 c 

188 c thermal effects option - temp from load curve 1 

189 c 

190 if (%isimjtnode.ne.5) then 

191 teo 1 

192 endif 

193 c 

194 tssfO.O 

195 c 

196 c print initial time step size 



198 cprtflg 1 FIG. 8E 

199 c 

200 c .... turn off (0) or on (1) SAND database flag 

201 c 

202 edsdfO 

203 c 

204 nrest 90000 

205 nrunr 95000 ; 

206 c 

207 c .... DYNA3D stent to compression cyl 

208 c 

209 sid 1 dni 

210 c sfif 

211 cmfif 

212 pnlts 1.0e-0 

213 pnltm 1.0e-0 

214 ; 

215 c 

216 c .... DYNA3D tied interface to interconnects if multisegment 

217 c 

218 if(%imodel.eq.0)then 

219 sid 2 tied 

220 ; 

221 endif 

222 c 

223 c .... end DYNA3D commands 

224 c 

225 endif 

226 c 

227 c *** NTKE3D Analysis Options *** 

228 c 

229 if (%inike.eq. 1 ) then 

230 echo Making NIKE3D input file . . . 

231 nike3d 

232 nikeopts 

233 c 

234 c .... temperatures follow load curve 1 

235 c ** manually add tref=l .0 on matl 2 control card cols 26-35 ** 

236 c 

237 teo 1 

238 c 

239 if (%isim_mode.eq.5) then 

240 anal dyn 

241 neig20 

242 shift 69 

243 iplt 1 

244 nsbrr 1 

245 stifcore 1 






246 
247 



bfgscore 
bwmo new 
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248 c 

249 c element constitutive data incore 

250 c 

251 bforlO 

252 sfor 10 

253 bef 11 

254 c 

255 c .... linear solver 

256 c 

257 lsolver fissle 

258 c 

259 elseif (%isim_mode.ne.5) then 

260 c 

261 c .... time step analysis 

262 c 

263 nstep 100 

264 delt 0.0100 

265 anal stat 

266 c 

267 c .... step tol of le-2 is OK for predel compression 

268 c 

269 if (%isim_mode.eq. 1 .or.%isim_mode.eq.2) then 

270 dctol -1.0e-3 

271 elseif (%isim_mode.eq.3) then 

272 dctol -1.0e-2 

273 endif 

274 c 

275 c .... max iterations per stiffness reform 

276 c 

277 nibsr 20 

278 c 

279 c .... max stiffness reforms per step 

280 c 

281 msrf20; 

282 if (%isim_mode.eq. 1 .or.%isim_mode.eq.2) then 

283 iprt 1 

284 elseif (%isim_mode.eq.3.or.%isim_mode.eq.4) then 

285 iprt 25 

286 endif 

287 iplt 1 

288 nsbrr 1 

289 stifcore 1 

290 bfgscore 

291 bwmo new 

292 echo Bandwidth minimization ACTIVATED with "NEW" option 

293 c 

294 c element constitutive data incore 




296 bfor 10 FIG. 8G 

297 sfor 10 

298 bef 11 

299 c 

300 c .... linear solver 

301 c 

302 lsolver fissle 

303 c 

304 c .... solid element stent contact surface 

305 c 

306 sid 1 sv 

307 c 

308 if (%isim_mode.eq. 1 ) then 

309 c 

310 c .... below changed for sharp-edge laser-cut stent 

311 c 

312 pnlt 1.0e-3 

313 elseif (%isim_mode.eq.2) then 

314 pnlt 0.01 

315 elseif (%isim_mode.eq.3) then 

316 c 

317 c . ... essential to cut penalty for laser-cut stent predel compression 

318 c 

319 pnlt 0.001 

320 elseif (%isim_mode.eq.4) then 

321 pnlt 1.0e-3 

322 c iaug 1 ; 

323 endif 

324 ; 

325 c 

326 c .... end block for time step only analysis 

327 c 

328 endif 

329 c 

330 c .... slidesurface between interconnects and segments 

331 c 

332 sid 2 tied 

333 ; 

334 c 

335 c .... slidesurface between stent and aortic wall 

336 c 

337 if (%isim_mode.eq.6) then 

338 echo *** Add activation time of 0.5 to slidesurface 2 *** 

339 sid 3 sv 

340 ; 

341 endif 

342 c 

343 c .... NIKE3D shell geometric stiffness (HL only) 





3V 



344 
345 



c 



segs 1 ; 



FIG. 8H 



346 c 

347 c . . . . end NIKE3D section 

348 c 

349 endif 

350 c 

351 c . . . . symmetry planes (omit for freq analysis) 

352 c 

353 if (%isim_mode.ne.5) then 

3 54 if (%isym.eq. 1 ) then 

355 c 

356 c .... Symmetric Model 

357 c 

358 c plane 1 

359 c 0.0 0.0 0.0 

360 c [-sin(60)] [-cos(60)] 0.0 

361 c 0.0005 symm ; 

362 c plane 2 

363 c 0.0 0.0 0.0 

364 c [-sin(60)] [cos(60)] 0.0 

365 c 0.0005 symm ; 

366 c 

367 else 

368 c 

369 c .... symmetry planes to remove rigid body modes for full model 

370 c 

371 plane 1 

372 0.0 0.0 0.0 

373 1.0 0.0 0.0 

374 .0005 symm ; 

375 plane 2 

376 0.0 0.0 0.0 

377 0.0 1.0 0.0 

378 .0005 symm ; 

379 endif 

380 endif 

381 c 

382 c 

3 83 if (%inike.eq.0) then 



. Load Curves for DYNA3D **** ADD DR FLAG TO INPUT FILE **** 



387 if (%isim_mode.eq. 1) then 

388 c 

389 c .... radial force 

390 c 

391 led 1 

392 0.O00E+00 1.000E+00 



c 



c 






393 
394 



7.500E-03 2.250E+02 
1.000E-00 2.250E+02 ; 



FIG. 81 



395 elseif (%isim_mode.eq.2) then 

396 c 

397 c .... flat plate compression, led 1 not used (dummy definition) 

398 c 

3 99 echo ! ! ! Flat plate not implemented for DYNA3D ! ! ! 

400 quit 

401 c 

402 elseif (%isim_mode.eq.3) then 

403 c 

404 c .... predelivery compression strain - 0.87 in. dia compressed to 12F 

405 c [check x-displ of stent center node to verify] 

406 c 

407 led 1 

408 O.O00E+00 1.000E+00 

409 1.000E-02 1.008E+03 

410 1.000E-00 1.008E-K)3 ; 

41 1 elseif (%isim_mode.eq.6) then 

412 c 

413 c .... anatomy deployment 

414 c (LC from radial comp) 

415 c 

416 led 1 

417 O.0O0E+O0 1.000E+00 

418 7.500E-04 1.000E+03 

419 9.000E-04 1.000E+03 

420 1.500E-03 1.000E+00 

421 1.000E-00 1.000E+00 ; 

422 endif 

423 c 

424 c .... load curve #2 only used for flat plate compression 

425 c 

426 led 2 

427 0.000E+00 0.000E+00 

428 1.000E+00 0.000e-00 ; 

429 endif 

430 c 

43 1 if (%inike.eq. 1 ) then 

432 c 

433 c .... ******** Load Curves for NIKE3D ********** 

434 c 

435 if (%isim_mode.eq. 1 ) then 

436 c 

437 c .... radial force 

438 c 

439 led 1 

440 O.OOOE+00 1.000E+00 

44 1 1 .000E+00 3 .000E+02 ; 



442 elseif (%isim_mode.eq.2) then TTT/^ O T 

443 c rlvjr. oJ 

444 c .... flat plate compression, led 1 not used (dummy definition) 

445 c 

446 led 1 

447 0.000E+00 1.000E+00 

448 1.000E+00 0.000E+00 ; 

449 elseif (%isim_mode.eq.3) then 

450 c 

45 1 c .... predelivery compression strain - 0.87 in. dia compressed to 12F 

452 c [check x-displ of stent center node to verify] 

453 c 

454 led 1 

455 0.000E+00 1.000E+00 

456 1.000E+00 1.008E+03 ; 

457 elseif (%isim_mode.eq.4) then 

458 c 

459 c .... initial expansion strain - 4/5 mm OD to 1 5/27 mm OD 

460 c [check x-displ of stent center node to verify] 

461 c 

462 led 1 

463 c .... thermal load (activate TEO above) 

464 O.OOOE+00 1.000E+00 

465 1.000E+00 -1.008E+03 ; 

466 c .... prescribed displacement 

467 c 0.000E+00 O.O00E+00 

468 c 1.000E+00 1.000E-01 ; 

469 c 

470 elseif (%isim_mode.eq.5) then 

471 c 

472 c .... must define load curve since TEO active even if unused for freq 

473 c 

474 c .... initial expansion strain - 4/5 mm OD to 1 5/27 mm OD 

475 c [check x-displ of stent center node to verify] 

476 c 

477 led 1 

478 c .... thermal load (activate TEO above) 

479 O.000E+00 1.000E+00 

480 1.000E+00-1.008E+03 ; 

48 1 elseif (%isim_mode.eq.6) then 

482 c 

483 c .... anatomy deployment - 0.87 in. dia compressed to 12F 

484 c 

485 led 1 

486 O.OOOE+00 1.000E+00 

487 0.500E+00 5.000E+02 

488 1.000E+00 1.000E+00 ; 

489 endif 

490 endif 



/ HL( 



491 c mp OT7 

492 c stent parts 1 • OIV 

493 c 

494 include irss.tg 

495 c 

496 c anatomy parts 

497 c 

498 if (%isim_mode.eq.6) then 

499 c 

500 c .... convert anatomy data from cm to inch units 

501 c 

502 control 

503 csca[L/2.54] 

504 c 

505 c .... import meshed anatomy data for stent deployment 

506 c (this is an aortic stent) 

507 c 

508 include tpeg.part_ct_aorta3 

509 csca 1.0 

510 merge 

511 if (%inike.eq.l) then 

512 c 

513 c .... set material properties for aortic wall 

514 c 

515 include aorta.materials_nike 

516 endif 

517 endif 

518 c 

519 c stent materials 

520 c 

521 if (%inike.eq.l) then 

522 if (%isim_mode.eq. 1 .or.%isim_mode.eq.2) then 

523 include istent.mats_nike_solid 

524 echo NiTi model for radial force/flat plate analysis 

525 elseif (%isim_mode.eq.3) then 

526 include istentmats_compress_nike_solid 

527 echo NiTi model for predelivery compression strain 

528 elseif (%isim_mode.eq.4) then 

529 include istent.mats__compress_nike_solid 

530 echo NiTi model for initial expansion strain 

53 1 elseif (%isim_mode.eq.5) then 

532 include istent.mats_nike_freq_solid 

533 echo NiTi model for frequency analysis 

534 elseif (%isim__mode.eq.6) then 

535 include istent .mats_nike__solid 

536 echo NiTi model for anatomy deployment 

537 endif 

538 c 

539 elseif (%inike.eq.0) then 



J// 



540 if (%isim_mode.eq. 1 ,or.%isim__mode.eq.2) then T?TS~^ O T 

541 include istent.mats_dyna_solid JtIvJ. Ol_/ 

542 echo NiTi model for radial force/flat plate analysis 

543 elseif (%isim_mode.eq.3) then 

544 include istentmats_compress_dyna_solid 

545 echo NiTi model for predelivery compression strain 

546 elseif (%isim_mode.eq.4) then 

547 include istentmats_compress_dyna_solid 

548 echo NiTi model for initial expansion strain 

549 elseif (%isim_mode.eq.6) then 

550 include istent.mats_compress_dyna_solid 

551 echo NiTi model for anatomy deployment 

552 endif 

553 endif 

554 c 

555 c .... cylindrical compression for radial force or predelivery compression 

556 c 

557 if (%isim_mode.eq. 1 .or.%isim_mode.eq.3 .or.%isim_mode.eq.4.or.%isim_mode.eq.6) then 

558 c 

559 if (%isym.eq. 1) then 

560 include cylinder.parts_sym 

561 else 

562 include cylinder.parts 

563 endif 

564 endif 

565 c 

566 if (%inike.eq.l) then 

567 include cylinder.materials_nike 

568 elseif (%inike.eq.O) then 

569 include cylinder.materials_dyna 

570 endif 

571 c 

572 stp .0001 

573 c 

574 c .... Constrain stent node(s) in z-direction for time-hist analysis 

575 c 

576 if (%isim_mode.ne,5) then 

577 merge 

578 c 

579 c .... nset for 3-segment model 

580 c nset zconstr = 1 8149 8687 9215 9747 ; 

581 c echo ** Bottom 12-crown node list Constrained in Z-translation ** 

582 c 

583 c .... nset for 6-crown only 

584 echo ** Bottom 6-crown node list constrained in z-dir ** 

585 nset zconstr - 1 43 97 151 448 ; 

586 b nset zconstr dz 1 ; 

587 endif 

588 c 
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